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Laser-Doppler velocimeter measurements were made in the flow field of a stalled 
cascade of controlled-diffusion stator blades. Tests were conducted at 10 degrees of 
incidence above the design inlet-flow angle in order to verify previous measurements. A 
unique method of data presentation was offered to characterize the unsteady positive and 
negative velocities about their mean. Laser-sheet flow visualization was performed and 
showed the unsteadiness of the stalling phenomena within the blade row passages. 
Additionally, a leading-edge separation bubble was observed and reversed flow was 
measured within the bubble using the laser-Doppler anemometer. Power-spectral density 
and auto-correlation analysis results are presented. In parallel with the experimental 
measurements, a computational fluid-dynamics study was initiated in an attempt ultimately 
to predict stall. Viscous flow computations were completed at on-design and near-stall 
conditions. Pressure distributions, separation bubble re-attachment locations, and loss 
predictions were compared to previous experimental studies. Generally good agreement 
was obtained for on-design conditions with degradation in prediction nearing stall. 



v 



TABLE OF CONTENTS 



I. INTRODUCTION 1 

II. EXPERIMENTAL SETUP 3 

A. LOW-SPEED CASCADE WIND TUNNEL 3 

B. INSTRUMENTATION 3 

C. EXPERIMENTAL PROCEDURE 6 

1 . Surface Pressure Measurement 6 

2. LDV Measurements 6 

a. Inlet Surveys 10 

b. Passage Surveys 10 

c. Wake Surveys 10 

d. Laser Power-Spectrum / Auto-Correlation Measurements 10 

3. Flow Visualization 1 1 

III. COMPUTATIONAL FLUID DYNAMICS 13 

A. GRID GENERATION 13 

B. COMPUTATIONAL FLOW SOLVER 15 

IV. RESULTS AND DISCUSSION 17 

A. BLADE-SURFACE PRESSURE COEFFICIENT 17 

B. LDV SURVEY COMPARISONS 18 

1. Inlet Survey (Station 1) 18 

2. Passage Survey (Station 7) 18 

3. Wake Survey (Station 19) 23 

C. PASSAGE SURVEY (STATION 3) 26 

D. POWER-SPECTRUM / AUTO-CORRELATION ANALYSIS 33 

E. FLOW VISUALIZATION 36 



vii 



F. COMPUTATIONAL RESULTS 



38 



1. Comparison of Blade-Surface Pressure Distributions........... 38 

2. Leading-Edge Separation Bubble 43 

3. Losses 46 

V. CONCLUSIONS AND RECOMMENDATIONS 48 

APPENDIX A: EXPERIMENTAL-MEASUREMENTS LOG 51 

APPENDIX B: CALCULATION OF INLET-FLOW REFERENCE VELOCITY 53 

APPENDIX C: VELOCITY HISTOGRAMS, TIME HISTORY PLOTS, POWER 
SPECTRUM PLOTS, AND AUTO-CORRELATION PLOTS 56 

1. Station 3, Point 7 56 

a. 3K Data Points 56 

b. 30K Data Points 60 

2. Station 7, Point 3 64 

a. 3K Data Points 64 

b. 30K Data Points 68 

c. Ganaim 3K Data Points 72 

3. Station 7, Point 20 76 

a. 3K Data Points 76 

b. 3 OK Data Points 80 

c. Ganaim 3K Data Points 84 

APPENDIX D: REDUCED DATA FROM STATIONS 1, 3, 7, 19 88 

APPENDIX E: CFD GRID GENERATION 92 

APPENDIX F: RVCQ3D CODE INPUT 94 

APPENDIX G: DENSITY RESIDUALS FOR CFD SOLUTIONS 95 

APPENDIX H: CP VERSUS X/C FORTRAN PLOTTING PROGRAM 96 

APPENDIX I: UPSTREAM FLOW ANGLE FORTRAN PROGRAM 99 

REFERENCES 101 

INITIAL DISTRIBUTION LIST 103 



vm 



LIST OF FIGURES 



1 . Schematic of the Cascade Wind Tunnel 4 

2. Laser-Sheet Flow- Visualization System 5 

3. LDV Survey Stations Between Blades 6 and 7 7 

4. LDV Power-Spectrum / Auto-Correlation Locations 9 

5. Two-Dimensional Periodic C-Grid 14 

6. Pressure Distribution Comparison 17 

7. Velocity Comparison at Station 1 19 

8. Turbulence Intensity Comparison at Station 1 19 

9. Velocity Comparison at Station 7 20 

10. Turbulence Intensity Comparison at Station 7 20 

11. U Velocity Comparison at Station 7 22 

12. V Velocity Comparison at Station 7 22 

13. Velocity Comparison at Station 19 24 

14. Turbulence Intensity Comparison at Station 19 24 

15. U Velocity Comparison at Station 19 25 

16. V Velocity Comparison at Station 19 25 

17. U Velocity Survey at Station 3 28 

18. U Velocity Histograms at Station 3 29 

19. V Velocity Survey at Station 3 30 

20. V Velocity Histograms at Station 3 31 



IX 



21. U-V Velocity Correlation Plot 32 

22. Power Spectrum for 3 OK Points 34 

23. Power Spectrum for 3 K Points 35 

24. 40 Degree Pressure Distribution 40 

25. 43 Degree Pressure Distribution 40 

26. 46 Degree Pressure Distribution 41 

27. 48 Degree Pressure Distribution 41 

28. Non-Dimensional Pressure Contours for 48 Degrees 42 

29. Leading-Edge Velocity Plot 44 

30. Separation-Bubble Reattachment 45 

3 1 . Loss Comparison 47 



x 



LIST OF TABLES 



1. LDV Survey Stations and Locations 6 

2. LDV Power-Spectrum / Auto-Correlation Locations 8 

3. Spectral- Analysis Maximum Frequencies 36 



xi 



Xll 



ACKNOWLEDGMENTS 



I would like to thank Dr. Garth Hobson for his tireless enthusiasm towards my 
work. It kept me motivated throughout. Thanks to Dr. Raymond Shreeve for first 
putting the curiosity there. And to Rick Still and Ted Best for their technical expertise and 
support. Finally, to Dawn, my wife and best friend, who supported what I believed to be 
important work while tending to far greater work at home ... our girls, Jessica and Sydney. 



xiii 



I. INTRODUCTION 



Compressor stall and off-design behavior impose performance constraints on gas 
turbine engines. The open literature currently lacks adequate documentation of 
quantitative measurements of turbomachinery flow fields at and beyond stall. This 
knowledge is essential for predicting the onset of stall with computational methods and 
optimizing blade geometries for off-design range with minimal degradation in design-point 
performance. Towards optimum design-point performance, controlled-difiusion (CD) 
blading was developed. The purpose of these blades was to avoid boundary layer 
separation and, in the case of supercritical (SC) blading, to be shock-free in the transonic 
range. 

Extensive laser-Doppler velocimeter (LDV) investigations [Refs. 1 - 3] of flow 
through a controlled-difiusion compressor cascade have been conducted in the region 
below stall at the Naval Postgraduate School (NPS). Recently, the reported cascade was 
stalled with an inlet-flow angle of 50 degrees and an extensive mapping of the flow field 
was performed using a two-component LDV system by Ganaim [Ref. 4], The motivation 
for the present study was to partially map the flow field near stall and to establish the 
repeatability of the results reported in Reference 4. Between Ganaim’s work and the 
present, the tunnel was reconfigured twice to pre-stall flow angles. This was done to 
confirm the repeatability of the flow field and to expand the available database for 
unsteady separated flows. Also, since the previous data were presented on a mean-flow 
basis, it was desired to more properly characterize the mean flow into its positive and 
negative components. This required a special method of post-processing the velocity- 
histogram data, as described within this report. 

It was also desired to characterize the unsteady flow using power-spectrum and 
auto-correlation analyses. It was hypothesized that there could exist a fundamental 
frequency for the stalled cascade, that it could be determined from the measurements made 
in the present study and should match that of the previous data [Ref. 4], 

In parallel with the LDV study, a computational fluid dynamics (CFD) study was 
initiated in an attempt ultimately to predict stall. Previous CFD studies on the present 
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cascade include those of Elazar [Ref. 1], Stow, et al [Ref 5], Hobson [Ref. 6], and Kang, 
et al [Ref. 7], 

Before attempting to calculate the unsteady, stalled, flow through a cascade (as 
was measured in the present study), it was necessary first to evaluate the code’s predictive 
ability at incidence angles off-design and near the stall point. A quasi-three-dimensional 
computer code [Refs. 8-9] for the solution of the thin-layer Navier-Stokes equations was 
used to predict the cascade performance from design (40 degrees inlet-flow angle) to near 
stall (48 degrees inlet-flow angle). The code was used to successfully compute the blade 
surface pressure distribution, the leading edge separation bubble and total pressure loss 
coefficient, for inlet flow angles from 40 to 48 degrees. 
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II. EXPERIMENTAL SETUP 



A. LOW-SPEED CASCADE WIND TUNNEL 

The subsonic cascade wind tunnel contained 20 controlled-diffusion blades. Each 
blade was 10 inches in span, 5.01 inches in chord and separated by 3 inches in the 
pitchwise direction. The cascade solidity (blade chord-to-pitch ratio) was 1.67 and 
Reynolds number based on chord was approximately 700,000. A detailed description of 
the test facility, test section and CD blading was fully documented in Reference 1 . Figure 
1 shows a schematic of the cascade. The tunnel was configured for an inlet-flow angle 
Oi) of 50 degrees, and an exit-flow angle (Pz) such that the downstream end-wall static 
pressure distribution was uniform. 

B. INSTRUMENTATION 

Blade surface pressure measurements were recorded using a 48-channel Scanivalve 
system controlled by an HP-9000 computer. Figure 1 shows the “fully-instrumented 
blade” location within the test section. Reference 1 gives the location of the blade surface 
pressure points. A full description of the data acquisition system is given in Reference 10. 

A two-component, four-beam TSI Model 9100-7 LDV system was used for all 
flow field measurements. The LDV setup, including laser type, optics, atomizer, seeding, 
traverse mechanism and data acquisition were thoroughly documented in Reference 4. 

The flow was visualized by injecting fog particles of 2 pm in size through a wand 
inserted through the end wall upstream of the blading. A 300 mW air-cooled argon-ion 
laser was used as the light source for laser-sheet flow visualization. The laser beam was 
transmitted to a lens system via a fiber-optic cable. The lens system ended with a 
cylindrical-piano lens, which produced the laser sheet. The laser sheet was introduced into 
the cascade tunnel, ahead of the blades, through the sidewall. A schematic of the system, 
showing the location of the introduction of the fog and the position of the video camera 
(located outside the window), are shown in Figure 2. 
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SIDEWALLS 




Figure 1. Schematic of the Cascade Wind Tunnel 
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LOCATION 




Figure 2. Laser-Sheet Flow- Visualization System 
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c. 



EXPERIMENTAL PROCEDURE 



1. Surface Pressure Measurement 

The tunnel plenum chamber pressure was set at 12 inches of water (gauge) and the 
blade surface pressure measurements were taken as described in Appendix B of Reference 
10. Inlet total and static pressures were recorded 2 1/3 chord lengths ahead of the blades. 
These pressures were used to non-dimensionalize the blade-surface pressure 
measurements as a coefficient of pressure; i.e., C p = (p loca i - p«) / (p t * - p®). 

2. LDV Measurements 

A full listing of the experimental runs is contained in Appendix A. Blade-to-blade 
surveys were conducted at eleven stations ahead of, between, and downstream of the 
blade row. All surveys were taken by collecting 3,000 data points at each location in the 
survey. All LDV data were collected in the “coincidence mode” with a 2.0x1 0 5 psec 
window used for validation. Table 1 summarizes the survey stations and locations. Figure 
3 illustrates the same information graphically. The numbers in parentheses on the figure 
indicate the horizontal distance between measurement points. The present study gave an 
average inlet-flow angle of 50.05 degrees. Reference 4 results were for an average Pi of 
50.21 degrees. 



Station 


Region 


Location (% Axial Chord) 


1 


Inlet 


30.9 % forward of LE (1.5 inches forward of LE) 


lb 


a 


4.28 % “ 


Id 


tt 


4.28 % “ 


2a 


Passage (Suction Side) 


0 % (a), LE 


3 


a 


5.15% aftofLE 


7 


a 


30.9% “ (1.5 inches aft of LE) 


14 


a 


93.6 % “ 


15 


a 


98.7 % “ 


17 


Wake 


6.18% aftofTE 


1 18 


It 


12.7% “ 


19 


a 


20.6% “ (1.0 inches aft of TE) 



Table 1. LDV Survey Stations and Locations 
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X (inches) 



Figure 3. LDV Survey Stations Between Blades 6 and 7 
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Additionally, eleven specific locations were measured, with an increased number of 
data points, for power-spectrum and auto-correlation analyses. Table 2 below summarizes 
the measurements and the number of data points gathered for each point. Figure 4 
depicts the stations at which these dedicated runs were performed. The points close to the 
blade were chosen to be within the intermittent reverse-flow region, and the outboard 
points were chosen to be outside the reverse-flow region. 



Station 


Point 


Location 


Number of Points Acquired 


1 


7 


(0.0000, -6.2920) 


3K, 30K, 50K ( 


2a 


1 


(-1.1331,-4.7905) 


3K, 30K 


3 


7 


(-0.8394, -4.5420) 


3K, 30K | 


7 


3 


(-0.1225, -3.2920) 


3K, 30K 


7 


20 


(0.3654, -3.2920) 


3K, 30K 


10 


15 


(0.3988, -1.7921) 


30K 


10 


25 


(1.0000, -1.7921) 


30K 


17 


23 


(0.2500, 0.3620) 


30K 


17 


31 


(1.5000, 0.3620) 


30K 


19 


23 


(0.2500, 1.0620) 


30K 


19 


31 


(1.2500, 1.0620) 


30K 



Table 2. LDV Power-Spectrum / Auto-Correlation Locations 



The tunnel reference velocity (V re f) for each run was computed as described by 
Elazar in Reference 1. The results are shown in Appendix B. The reference velocity was 
used to non-dimensionalize the data for comparison with previous studies. 

Prior to each run, the probe volume was always referenced to the same location 
between the blades using an alignment tool. This allowed for a high confidence level in 
specifying station position for comparison with previous studies. The alignment procedure 
was developed and described in detail in Reference 1 . 
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X (inches) 



Figure 4. LDV Power-Spectrum / Auto-Correlation Locations 
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a. Inlet Surveys 

The inlet region was surveyed across stations 1, lb, and Id (Figure 3). The 
survey spanned three and one-third passage widths in the pitchwise direction (-9 inches to 
+1 inch) for station 1, and one and one-third passage widths in the pitchwise direction (-2 
inches to + 2 inches) for stations lb and Id. The laser beam system was horizontal when 
measuring at station 1 . A 4 degree pitch-up of the LDV table was required for stations 
lb and Id to avoid blade shadowing of the vertical beam, and to achieve closer access to 
the leading edge. The potential problems due to pitching and yawing the LDV system 
were treated by Hobson and Shreeve in Reference 3 . Their analysis showed the maximum 
spatial error due to probe volume orientation to be 0.3 mm, which is the minimum 
diameter of the measurement volume. 

The optics were configured the same throughout the present study with the 
488-nm blue beam measuring the horizontal (U) component and the 514.5-nm green beam 
measuring the vertical (V) component. Frequency shifting was performed as outlined in 
Reference 4 in order to capture the large reversed-flow velocities. 

b. Passage Surveys 

The passage surveys were performed between blades 6 and 7. All passage 
surveys extended 1.8 inches (approximately 60 % passage width) away from the suction- 
side surface of blade 6. The LDV system was pitched up 2 degrees and yawed 4 degrees 
left for stations 2a, 3 and 7 for better access to the suction side of the blade. Stations 10, 
14 and 15” were measured with only 4 degrees of left yaw. 

c. Wake Surveys 

The LDV system was reset to zero pitch and zero yaw for measurements at 
stations 17, 18, and 19. The wake surveys were centered around blades 6 and 7 and 
spanned two passage widths (6 inches). 

(L Laser Power-Spectrum /Auto-Correlation Measurements 

It was desired to expand the analysis of several points throughout the flow 
field by collecting a larger amount of data points. Instead of the standard 3,000 points 
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gathered at each location, 30,000 data points were chosen for the locations listed in Table 
2. This choice corresponded to the limiting number of data points that could be acquired 
at positions with low data rates such that the acquisition software reached its maximum 
timeout of 30 minutes. This corresponded to a minimum data rate of 16.7 Hz in order to 
collect 30,000 samples. Stations 2a and 3, near the separation-bubble region, were most 
affected by this limitation. 

3. Flow Visualization 

The laser-sheet flow visualization was photographed through the optical access 
window using a hand-held 8mm camera at 30 frames/second and a VHS camera at 60 
ffames/second. The filming was performed at night to reduce glare and optimize the flow 
field image contrast. A total of approximately 15 minutes of footage was obtained of the 
intermittent reversed flow inside the blade passage and of the flow near to the passage 
entrance, including the leading-edge separation bubble. 
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III. COMPUTATIONAL FLUID DYNAMICS 



A. GRID GENERATION 

The grid generation program employed was a modified version of GRAPE (GRids 
about Aerofoils using Poisson’s Equation), [Ref. 11]. This version allowed for the 
generation of a periodic C-type grid for turbomachinery cascades. The namelist input file, 
along with a description of the effort put forth in refining the leading edge geometry for 
the grid, is described in Appendix E. 

The grid, a portion of which is depicted in Figure 5, contained 240 grid points 
wrapped around the blade in the streamwise direction and 49 grid points in the blade-to- 
blade direction. The grid was much finer in the normal direction near the surface for 
improved viscous flow resolution, particularly of the leading edge separation bubble. The 
outer, or periodic, boundaries corresponded in size to the cascade wind tunnel as 
described in Experimental Setup. The inlet region was located 30% of an axial chord 
upstream of the leading edge. The wake region extended two chord lengths downstream 
of the trailing edge. 
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Figure 5. Two-Dimensional Periodic C-Grid 
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B. COMPUTATIONAL FLOW SOLVER 



Numerical computations were made using Rotor Viscous Code Quasi-3D 
(RVCQ3D) of Reference 8. RVCQ3D solved the thin-layer Navier-Stokes equations by 
employing an explicit, multi-stage, Runge-Kutta algorithm. It could account for the 
effects of rotation, radius change and stream surface thickness variation in blade-to-blade 
turbomachinery flows. The Baldwin-Lomax algebraic eddy- viscosity model [Ref. 12] was 
used to characterize the turbulence within the flow field. 

Computations were performed at the following inlet flow angles: 40, 43, 46, and 
48 degrees. This was done to allow comparisons with existing experimental results. 
Validation of the computed results was achieved by comparing pressure coefficient 
distributions around the blade surface, leading-edge separation bubble sizes, and loss 
computations. 

A typical namelist input file is shown in Appendix F. All solutions were run to 
6,000 iterations and axial-velocity-density ratio (AVDR) was held at 6% for all 
computations. The AVDR was modeled by defining two stream surface points at the inlet 
and exit using the namelist 6 parameter ‘NMN’ as shown in Appendix F. 

Although the inlet flow angle was specified initially, the flow solver computed to a 
new inlet condition based on the prescribed exit-static to inlet-stagnation pressure ratio. 
The inlet flow angle varied due to changes, primarily, in the tangential velocity 
component, as the axial velocity did not change due to continuity constraints. 

Solutions were computed using the National Aerodynamic Simulation (NAS) 
Program CRAY C90 supercomputer and the NPS CRAY Y-MP EL98 supercomputer. 
CPU run time for the NAS C90 was between 3-4 minutes and the NPS EL98 required 
about 55 minutes. Analysis of the solution files was carried out using FAST graphics 
software [Ref. 13] and separate FORTRAN programs given in Appendices H (blade- 
surface pressure distribution) and I (computed inlet-region flow angle). 
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IV. RESULTS AND DISCUSSION 



A. BLADE-SURFACE PRESSURE COEFFICIENT 

A comparison of pressure coefficient around the blade surface between the present 
study and Reference 4 is shown in Figure 6. The results, after reconfiguring the tunnel, 
showed reasonable agreement with the pressure distributions obtained earlier at the 
leading edge; particularly on the pressure side. There was some difference near the 
trailing edge which was possibly due to a difference in setting of the tunnel’s exit-plane 
tailboards. The setting of the tailboards significantly affected the exit static-pressure 
distribution. 




Figure 6. Pressure Distribution Comparison 
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B. 



LDV SURVEY COMPARISONS 



The data analysis and comparisons which follow were performed without any 
refinement bounds for the velocity histograms. All the results from Reference 4 were re- 
analyzed in this manner. 

1. Inlet Survey (Station 1) 

Figures 7 and 8 show the non-dimensionalized mean velocity and turbulent 
intensities, respectively. The two studies showed variations in the mean velocity of less 
than 2% across three passage widths. The period in the variation in the mean velocity 
corresponded to the 3 inch blade separation, indicating that the blades were felt 30% of an 
axial chord ahead of the leading edge. The U and V turbulence-intensity levels were 
similar, and this indicated the turbulence was isotropic at the inlet. The increase in 
turbulence intensity from 2% to 4% over the extent of the survey is noted. This behavior 
was not observed at lower inlet-flow angles. 

2. Passage Survey (Station 7) 

Mean-flow results for station 7 are presented in Figure 9. There was good 
agreement across the passage which further established repeatability with Reference 4’s 
results. The turbulence intensity in Figure 10 also agreed well and showed anisotropic 
behavior of the turbulence near the blade surface out to 50 % pitch. Beyond this point the 
turbulence became isotropic in the inviscid region, as the levels of turbulence in the U and 
V components were comparable. The data were presented as turbulence intensities, yet 
because the flow field was intermittently reversible the standard deviation of the 
histograms, given by the LDV, did not represent the turbulence in the flow. Rather, these 
should be looked upon as a measure of the unsteadiness in the flow, which was a 
combination of the random turbulence and intermittent reverse flow. 
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Station 1 Mean Velocity 
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Figure 7. Velocity Comparison at Station 1 



Station 1 Turbulence Intensity 
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Figure 8. Turbulence Intensity Comparison at Station I 
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Station 7 Mean Velocity 
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Figure 9. Velocity Comparison at Station 7 



Station 7 Turbulence Intensity 

AC\ 




oc 


















75 


OO 

on 


















ou - 
oc 


















S? 

• — ■ on 


















3 

15 - 


















in - 


















c 


















o 

n 


















u 

-0. 


25 0.1 


do o.; 


25 0. 1 . 


50 0. 

X (in 


75 1.1 

ches) 


do i.: 


25 U 


so 


— Ganaim U-Tu — b— Ganaim V-Tu Present Study U-Tu —■—Present Study V-Tu 



Figure 10. Turbulence Intensity Comparison at Station 7 
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Figures 1 1 and 1 2 separate the mean horizontal (U) and mean vertical (V) velocity 
components into their respective positive and negative mean values. The velocity 
histograms for each point across the passage were specially edited such that only negative 
velocities were retained and averaged to create a negative velocity trace and, likewise, 
only positive velocities were retained and averaged at each point to create a positive 
velocity trace. It was apparent from Figures 1 1 and 12 that the weighted average of the 
negative trace and positive trace would make up the mean velocity for each component 
across the passage as previously shown in Figure 9. 

This analysis showed more of the character of the intermittent reversed flow than 
did a plot of the mean velocity. What occurred in the unsteady stalling process was not 
represented by the arithmetic mean. The envelope about the mean contained the positive 
and negative velocities that truly characterized the flow’s behavior. Since the data were 
taken in coincidence mode, one disadvantage of this form of presentation was that during 
editing of the positive U velocity-component data, all of its negative points were 
discarded, as well as the corresponding V velocity-component data. The opposite but 
equivalent process occurred during editing of the V velocity component; such that all the 
data that were left were the velocity components in the first quadrant of a flow that was 
oscillating with components in all four quadrants. It was apparent from the figures that by 
30 % pitch, the flow no longer exhibited reverse velocities in either the U or V 
components. 
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Figure 11. U Velocity Comparison at Station 7 
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Figure 12. V Velocity Comparison at Station 7 
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3. Wake Survey (Station 19) 

The wake surveys were compared at station 19. Figures 13 and 14 show the 
results and indicate good general agreement for the mean velocity and turbulence 
intensity/unsteadiness across two passage widths for the two studies. The periodic 
influence of the blades was strongly evident in the V component of Figure 13. Also, the 
asymmetry of the wake was especially apparent with the larger velocity deficit on the 
suction side of the blade. Figure 14 once again points out the anisotropic behavior of 
turbulence/unsteadiness in this same region. 

Figures 15 and 16 are the results of the special editing as outlined for station 7 in 
Figures 11 and 12. It was of interest to note that the V-component reversed-flow 
velocities at station 19 were of the same magnitude as the V-component reversed-flow 
velocities at station 7. 
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Figure 13. Velocity Comparison at Station 19 
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Figure 14. Turbulence Intensity Comparison at Station 19 
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Station 19 U Velocity 
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Figure 15. U Velocity Comparison at Station 19 
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Figure 16. V Velocity Comparison at Station 19 
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C. PASSAGE SURVEY (STATION 3) 



Measurements obtained at 32 locations at station 3 are shown in Figures 17-20. A 
detailed analysis to resolve the unsteady characteristics of the flow was made at the three 
locations (1, 7 and 20) shown in Figure 17 (U component of velocity) and Figure 19 (V 
component of velocity). The histograms of the U and V components of velocity at the 
three locations are shown in Figure 1 8 and Figure 20, respectively. Point 1 was located 
just off the blade surface. This was inside the leading-edge separation bubble, evidenced 
by the negative mean velocity of all the data points in the histogram (Figures 18a and 20a). 
The second location (point 7) was at 2.5 % pitch. This point was within the intermittent 
reversed-flow region, as indicated by the bi-modal velocity histograms (Figures 18b and 
20b). The third location (point 20) was at 17 % pitch, and was outside the region of 
intermittent reversed-flow. 

Figure 1 7 shows the mean pitchwise velocity along with the separate averages of 
the positive and negative components. The progression from a negative mean velocity to 
a positive mean velocity was apparent with the crossover distinctly evident in the bi-modal 
histogram at point 7. It was of interest to note that the arithmetic mean for the flow rarely 
occurred. The negative and positive velocity peaks were predominant, which was 
indicative of the unsteady stalled flow. 

The streamwise and pitchwise velocity component behavior at station 3 appear to 
be similar because the blade surface was oriented at approximately 45 degrees to the two 
pairs of laser beams (X-Y axes for the velocity measurements). Figure 21 provides an 
alternate view of the previous data by plotting the pitchwise (U) and axial (V) velocity 
components together. 
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Figure 17. U Velocity Survey at Station 3 
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Figure 18. U Velocity Histograms at Station 3 
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Station 3 V Velocity 
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Figure 19. V Velocity Survey at Station 3 
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Figure 20. V Velocity Histograms at Station 3 
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Station 3, Point 1 
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Figure 21. U - V Velocity Correlation Plot 
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D. POWER-SPECTRUM / AUTO CORRELATION ANALYSIS 



Power spectral density and auto-correlation analyses were performed on the data 
collected as previously indicated in Table 2. Only station 3 (point 7) and station 7 (points 
3 and 20) have been included in Appendix C. Therein are contained the velocity 
histograms and time-history plots for these points. This analysis was performed using 
TSI’s FIND software, version 4.03. The settings chosen for the analysis were the 
following: Window Type: Hanning, Fast Fourier Transform (FFT) Method: Random 
Correlation, and Correlation Type: Random Correlation. In all cases the default lagged- 
time slot (Ax) was used in order to satisfy the Nyquist sampling criterion. 

It was initially believed that a large data set was required to properly carry out this 
type of analysis. The data sets acquired for station 3, point 7, included 3K and 3 OK data 
points. The comparison of the two sets showed similarity in their behavior. The same can 
be said for station 7, points 3 and 20, which were compared between 3K, 30K, and 
Ganaim’s 3K data points. The only point of recurring significance throughout the power 
spectrum and auto-correlation plots was in the region between 20-25 Hz as shown in 
Figure 22 (30K points, present study), and Figure 23 (3K points, Reference 4) for station 
7, point 3. No further interpretation of the results obtained was attempted due to the many 
limitations of discretely-sampled, randomly-arriving data points. Additionally, since this 
routine involved averaging the data over the lagged-time slot (Ax), there was some 
question whether the true character of the flow had been captured. 
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Figure 22. Power Spectrum for 30K Points 
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Figure 23. Power Spectrum for 3K Points 
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The use of larger data sets did not appear to be of significance in the final analysis. 
The primary limitation turned out to be the data sampling rate. The maximum frequency 
for spectral analysis was limited to 1/2 the data rate due to the Nyquist criterion, which 
prevents aliasing of the data. Table 3 shows the data rates and maximum frequencies for 
the spectral analysis plots of Appendix C. 



Location 


Number of Data Points 


Data Rate 
(Hz) 


Maximum 

Frequency 

(Hz) 


Station 3, Point 7 


3K 


39 


19 


« 


30K 


49 


24 


Station 7, Point 3 


3K 


82 


41 


w 


30K 


127 


63 


«c 


3K (Ganaim) 


606 


303 


Station 7, Point 20 


3K 


69 


34 


k 


30K 


71 


35 




3K (Ganaim) 


613 


306 



Table 3. Spectral-Analysis Maximum Frequencies 



E. FLOW VISUALIZATION 

Visualization of reversed-flow velocities within a blade passage was documented 
on VHS videotape. The tunnel was run at a plenum chamber pressure of 12 inches H 2 O to 
match the conditions during the LDV measurements. The shutter speeds for the 8mm 
recording- were varied from l/250th to l/1000th and l/4000th of a second with the rate of 
30 frames/second remaining constant. It was found that the slower shutter speed of 
1 /250th of a second was better able to capture the streaklines of the intermittent reversed- 
flow due to longer particle traces for each image. 

Upon replaying the video? at normal and slower speeds, several observations were 
made for both the present study and the flow visualization performed in Reference 4. The 
reversed-flow region appeared to be confined to the blade passage and did not extend 
upstream of the leading edge. The leading-edge separation bubble trapped and 
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accumulated seed particles, due to injection from the reversed flow. This seed within the 
bubble was able to sustain itself even after the fog was turned off. The bubble also 
appeared to pulsate in size in accordance with the intermittent reversed-flow. There was 
some three-dimensionality to the reversed flow evidenced by small spanwise velocity 
traces. Since the laser sheet was only projected between two blades, it was not possible to 
determine whether the cascade was experiencing intermittent stall from blade passage to 
blade passage, or if the entire cascade was stalling simultaneously. 
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F. 



COMPUTATIONAL RESULTS 



1. Comparison of Blade-Surface Pressure Distributions 

The predicted and experimental blade-surface pressure distributions are shown for 
40, 43, 46, and 48 degrees in Figures 24-27 respectively. Experimental data for the 40 
and 43 degree case were from Reference 14. The 46 degree case used Reference 15 
experimental data and the 48 degree case used Reference 3 data. The AVDR parameter 
was held constant at 6% for all solutions, which showed good agreement for the first three 
cases. The 48 degree prediction over-estimated the experimental values aft of 60% chord. 
This was most probably due to under-estimating the stream-tube thickness variation at this 
inlet flow angle. Increasing the AVDR, as inlet flow angle increased, improved the 
agreement of the computational results; however, this was not the case experimentally, as 
documented in References 14, 16, and 17. 

The computed static-pressure recovery far downstream of the trailing edge and 
that measured (from Reference 18) showed good agreement for the 40 and 43 degree 
cases. The computed value for 40 degrees was 0.313 as compared with a measured range 
of 0.32 - 0.34. The computed value for 43 degrees was 0.357 as compared with a 
measured range of 0.37 - 0.39. The computed static-pressure rise far downstream is 
relevant to accuracy in predicting the losses using viscous codes, and it was shown here to 
agree well with experiments near design incidence. 

For Pi = 48 degrees, the predicted flow field pressure contours between the blade 
rows are shown in Figure 28. The pressure field was continuous at the grid boundaries 
showing good periodicity; however, there was some indication of grid contamination of 
the pressure field at the trailing edge cut. The suction peak and stagnation point were 
both clearly evident. The static-pressure rise, or adverse-pressure gradient, towards the 
trailing edge on the suction side was also apparent. Investigation of the velocity field in 
this region showed that there was flow separation aft of mid-chord, followed by re- 
attachment prior to the blunt trailing edge. Similar behavior occurred for the 46 degree 
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solution. Experimentally this was not the case, as the flow remained attached on the aft 
section of the suction surface. This indicated that the code may be predicting incipient 
stall at 48 degrees (or at a lower inlet-flow angle of 46 degrees); whereas, experimentally 
this was only found at 50 degrees. It appeared that an unsteady solution may have been 
computed as evidenced by the climbing density residuals, which are shown in Appendix G, 
and the coalescence of the pressure contour lines throughout the flow field (Figure 28). 
The 40 and 43 degree solutions did not show any suction-side separation on the aft 
portion of the blade, in agreement with experiments. 
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Figure 24. 40 Degree Pressure Distribution 
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Figure 27. 48 Degree Pressure Distribution 
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48 Degrees 




Figure 28. Non-Dimensional Pressure Contours for 48 Degrees 
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2. Leading-Edge Separation Bubble 

A leading-edge separation bubble was predicted for all computational solutions. 
The location of the computed separation point was insensitive to varying inlet flow angle. 
The re-attachment point for each case was determined graphically. Figure 29 shows the 
bubble region for Pi = 48 degrees in which the recirculating flow was evident. A 
comparison of computed and experimental re-attachment points is shown in Figure 30. 
The experimental data was from Reference 19, obtained using the china-clay technique. 
The plot showed an under-prediction of greater than 50% for the re-attachment locations, 
which is typical of bubble computations. The Bald win-Lo max turbulence model does not 
account for freestream turbulence, which, as pointed out in Reference 3, is necessary in 
order to properly characterize the development of the blade surface boundary layer, or the 
occurrence of separation. The two most likely reasons for the shortcoming of the 
computations are firstly, the limitations of the Baldwin-Lomax turbulence model and 
secondly, the lack of grid refinement in the streamwise direction aft of the leading edge. 
The trend of increasing downstream distance for re-attachment, as inlet flow angle 
increased, appeared to be properly computed. 
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Figure 29. Leading-Edge Velocity Plot 
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3. 



Losses 



The final assessment of RVCQ3D was its loss-prediction capability. Figure 31 
plots measured losses as reported in References 14, 16 and 17 with those computed in the 
present study. The first two computed cases of 40 and 43 degrees were within 5% of the 
measured results. The 46 degree case showed a sharp departure from the measured 
results, falling short by 45%. An attempt to increase the AVDR parameter for this case 
led to a small increase in the loss prediction but with a marked degradation of the blade 
pressure distribution. 

The output file for each solution contained a listing of the Baldwin-Lomax 
turbulence-model behavior by printing values of turbulent eddy viscosity for each grid 
point. Upon interrogation of this output file, the behavior of the Baldwin-Lomax 
turbulence model was seen to exhibit certain characteristics. Firstly, for all cases the 
turbulence model turned on inside the separation bubble region. Secondly, for the 40 
degree case, the model did not turn on until two grid points inside the bubble region, 
indicating a laminar region in the forward section of the bubble. The model then turned 
off again at the re-attachment point and turned back on further downstream. Thirdly, for 
48 degrees the model turned on at the separation point and remained on past the re- 
attachment point to the trailing edge. The transition criteria used for all the computations 
was that once the turbulence model computed a turbulent eddy viscosity (px) which was 
14 times larger than the laminar viscosity (pl), this value of p T was used to compute an 
effective viscosity (peff = pl + Pt ). Otherwise the effective viscosity was simply the 
laminar viscosity. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



Experimental flow field measurements near stall were successfully repeated for a 
controlled-diffusion compressor cascade using a two-dimensional laser-Doppler 
velocimeter system. The nominal inlet flow angle achieved was 50.05 degrees at a 
constant Mach number of 0.25 and Reynolds number of 700,000. The measurements 
were useful in that they showed repeatability with previous work by Ganaim, and 
validated a test case for two-dimensional unsteady separated flows. 

Analysis of the data was presented through a unique editing process which 
captured the character of the intermittent reversed flow. This process more properly 
represented the stalling region as a combination of primarily positive and negative flow, 
since the time-mean velocity rarely occurred in the intermittent region. 

Power-spectrum and auto-correlation calculations were completed on the data for 
eleven locations. Interpretation of these results was limited to their general behavior due 
to averaging of the data by the software package. Three locations were chosen for 
analysis and the results were found to depend greatly on the data acquisition rate. When 
expecting to perform power-spectrum analysis on experimental data, it is important to 
ensure that the data rate is high by repeated cleaning and proper re-fitting of the optical- 
access window. 

Numerical computations showed good predictability of the blade pressure 
distribution up to 46 degrees (approximately 6 degrees above design incidence). The 48- 
degree case began to deviate by up to 50% from measured values aft of two-thirds chord. 
A case-by-case basis for adjusting AVDR may be necessary to bring the computed results 
into better agreement with the measurements at higher inlet flow angles. 

A leading-edge separation bubble was formed in all computations in accordance 
with experiments. The extent of the computed bubble fell short of that seen in 
measurements. Computed losses agreed well at low incidence angles. As the inlet flow 
angle increased above 46 degrees, predicted losses deviated by 45% from the 
measurements. The use of a constant AVDR for numerical solutions greatly simplified the 
process. This was found to be acceptable for analyzing near on-design inlet flow angles. 
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Future computations of the flow through the cascade with a quasi-3d code should 
focus on achieving the following: 

1 . Continue numerical solutions to determine the optimum AVDR as inlet-flow 
angle varies, and compare with those measured in the cascade. 

2. Quantify what effects adjusting AVDR to match pressure distributions has on 
the computed losses. 

3. Increase the grid size in both the I and J directions for better characterization of 
the flow near the blade surface and around the leading edge separation bubble. 

4. Replace the Baldwin-Lomax turbulence model with a two-equation turbulence 
model to account for ffeestream turbulence intensity. 

5. Perform time-accurate solutions at higher incidence angles to see if 
unsteadiness is being properly modeled. 
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APPENDIX A: EXPERIMENTAL-MEASUREMENTS LOG 



A detailed log of the atmospheric conditions, tunnel plenum chamber 
conditions and LDV traverse-mechanism geometry for each experimental run. 
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Plenum 
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(inH 2 0) 


Plenum 

Temp 

<*F) 


Prandtl 
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(in H 2 0) 


Atmos. 

Pressure 
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LDV 

Pitch/ 

Yaw 
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11 14s1 a 


1-41 


11/14/94 


17:00 


11.9 


65 


17.0 


14.70 


0/0 










17:30 
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66 


16.8 


- 






stl 3k 
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11/30/94 


7:40 


12.0 


58 


17.0 


14.80 


0/0 










7:50 
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59 
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- 






stl 30k 
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11/30/94 


7:55 


12.0 


59 


17.0 


14.80 


0/0 
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12.1 


- 


- 


u 






stl 50k 
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11/30/94 


8:45 


12.0 


60 


17.0 


14.80 


0/0 










8:55 


- 


- 


■ 


- 




1b 


1121s 1b 


1-17 


11/21/94 


16:35 


12.0 


61 


16.9 


14.79 


+4/0 










16:51 


12.1 


63 


u 


- 




Id 


1122s Id 


1-33 


11/21/94 


16:55 


12.1 


64 


16.9 


14.79 


+4/0 










17:19 


• 


65 


17.1 


• 




2a 


1122s 2a 


1-32 


1 1/22/94 


8:45 


12.0 


61 


16.8 


14.85 


+2/4L 










9:35 


11.9 


64 


16.6 


- 






st2a 3k 


1 


11/30/94 


10:28 


12.0 


64 


16.8 


14.80 


+2/4L 










10:35 


- 


- 


■ 


- 
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1 


12/5/94 


7:55 


12.1 


58 


17.1 


14.70 


+2/4L 










8:27 


12.2 


60 
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14 69 




3 
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7:50 


11.9 


60 
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+2/4L 










8:45 
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66 


16.7 


14.84 
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61 


16.8 
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+2/4L 










10:13 


- 


62 


- 


- 






st37 30k 


7 
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8:35 


12.0 


61 


17.1 


14.69 


+2/4L 










9:07 


- 


62 


- 


14.70 




7 


1123s 7 


1-32 


11/23/94 


9:00 


12.0 


67 


16.8 


14.84 


+2/4L 










9:24 


■ 


68 


“ 


■ 
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9:40 


12.1 


60 


17.0 
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9:45 


■ 


- 


- 


■ 
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3 
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63 


16.9 
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17.0 


• 






st7 3k 
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8:05 
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“ 


■ 


■ 


“ 






st7 30k 


20 
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8:12 


11.9 


58 
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14.75 


0/0 










8:25 


- 


60 


“ 


■ 




10 


stIO eda 


15 


12/6/94 


9:28 


12.1 


63 


17.0 


14.76 


0/0 










9:43 


- 


“ 


- 


- 






stIO out 
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8:53 


12.0 


62 


168 


14.75 


0/0 










8:58 


12.1 


■ 


■ 


• 




14 


1127s 14 


1-32 


11/27/94 


15:33 


12.0 


58 


17.1 


14.82 


0/4L 










16:00 


12.2 


62 


17.3 


■ 




15 


1128s 15 


1-32 


11/28/94 


8:00 


12.0 


60 


17.0 


14.84 


0/4L 










8:35 


12.1 


62 


- 


- 




17 


1128s 17 


1-49 


11/28/94 


8:52 


12.0 


62 


17.0 


14.84 


0/0 










9:48 


12.2 


64 


- 


- 






stl 7 in 


31 


12/6/94 


8:30 


11.9 


61 


16.9 


14.75 


0/0 










8:37 


- 


- 


- 


■ 






stl 7 out 


23 


12/6/94 


8:40 


11.9 


61 


16.8 


14.75 


0/0 










8:48 


■ 


62 


■ 


• 
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Station 


Filename 


Points) 


Date 


Time 


Plenum 
Pressure 
(in H20) 


Plenum 

Temp 

(oF) 


Prandtl 
Probe 
(in H20) 


Atmos. 

Pressure 

(P*i) 


LDV 

Pitch/ 

Yaw 


18 


1128s 18 


1-49 


11/28/94 


10:00 


11.9 


64 


17.0 


14.85 


0/0 










10:43 


“ 


66 




44 




19 


1128s 19 


1-49 


11/28/94 


10:47 


11.9 


66 


17.0 


14.84 


0/0 










11:40 


12.1 


68 


44 


44 






stl9 in 


31 


12/6/94 


9:12 


12.2 


62 


17.3 


14.76 


0/0 










9:24 


12.3 




44 


44 






stl9 out 


23 


12/6/94 


9:02 


12.1 


62 


17.1 


14.75 


0/0 










9:09 


12.1 


44 


17.2 


14.76 
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APPENDIX B: CALCULATION OF INLET-FLOW REFERENCE VELOCITY 



cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 



OUTPUT FROM PROGRAM CALIBRATE 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

LEAST SQUARES STRAIGHT LINE CURVE FIT IS USED 
TO DETERMINE TUNNEL CHARACTERISTICS AT DIFFERENT SPEEDS 

NEWTON S METHOD IS USED TO DETERMINE THE REFERENCE VELOCITY 
FROM THE RECORDED AMBIENT PRESSURE AND TUNNEL PLENUM 
PRESSURE AND TEMPERATURE 

BEGIN DETERMINING TUNNEL CHARACTERISTICS 
FROM THE FOLLOWING MEASURED VALUES 

AXIAL VEL. TANGENTIAL VEL AMBIENT PRESS. PLENUM PRESS. PLENUM TEMP. 

M PER SEC. M PER SEC. INCHES MERCURY INCHES WATER DEG. C. 



18.3100 


22.1350 


29.9702 


1.7000 


13.8888 


32.4070 


38.5280 


29.9702 


4.7000 


13.8888 


42.5440 


50.1820 


29.9702 


8.1000 


15.5555 


47.2920 


56.3530 


29.9702 


10.0000 


16.6666 


51.7870 


61.9370 


29.9702 


12.1000 


17.2222 


55.3660 


66.7010 


29.9702 


14.0000 


18.3333 



CALCULATED VALUES FOR THE TUNNEL CONFIGURATION 



TOTAL VELOCITY MACH NUMBER 



MACH NUMBER FUNCT. PRESSURE RATIO 



0.2872654391 0E+02 
0. 5034501 3984E+02 
0.65789247298E+02 
0. 7356761 4295E+02 
0.80734660079E+02 
0. 86685 739064 E+ 02 



0.3781 8777200E-01 
0.66279705382E-01 
0.86361 830486E-01 
0.963872351 72E-01 
0.1 056761 5343E+00 
0.11 3249251 18E+00 



0. 4988029497 4 E-02 
0.1 52071 92279E-01 
0.25620261 31 9E-01 
0.31 7667591 19E-01 
0. 3800396621 2E-01 
0.43463396462E-01 



0. 4344 3528269 E-02 
0. 1 1919477 085E-01 
0.20366465634E-01 
0.25024235973E-01 
0.30121 036781 E-01 
0. 346867 26827 E-01 



CALLING LEAST SQUARES SUBROUTINE 
TO DETERMINE THE PRESSURE RATIO AS A FUNCTION OF MACH NO. PARAM 



PRESSURE RATIO = A1 * ANUX + A0 

MATRIX EQUATION 

0.60E+01 0.16E+00 A0 0.13E+00 

0.16E+00 0.53E-02 A1 0.42E-02 

A1 = 0.78899597626E+00 A0 = 0.16213280706E-03 
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REFERENCE CONDITIONS FOR EACH RUN 



AMB PRESS 
INCHES Hg 

29.9295 

30.1331 

30.1331 

30.1331 

30.1127 

30.1127 

30.2349 

30.1331 

29.9295 

30.1942 

30.1331 

29.9295 

30.2145 

30.1331 

29.9295 

30.1331 

30.0313 

30.0517 

30.0313 

30.1738 

30.2145 

30.2145 

30.0313 

30.0313 

30.2349 

30.2145 

30.0517 

30.0517 



PLENUM PRESS 
INCHES H20 

11.9000 

12.0000 

12.0500 

12.0000 

12.0500 

12.1000 

11.9500 

12.0000 

12.1500 

11.8500 

12.0000 

12.0000 

12.0000 

12.1000 

12.0000 

12.1000 

11.9000 

12.1000 

12.0500 

12.1000 

12.0500 

12.1000 

11.9000 

11.9000 

11.9000 

12.0000 

12.2500 

12.1000 



PLENUM TEMP 
DEGREES C 

18.3333 

14.4444 

14.4444 

15.5555 

16.6666 

18.3333 

16.9444 

17.7777 

14.4444 

17.2222 

16.6666 

16.6666 

20.0000 

15.5555 

17.2222 

14.4444 

14.4444 
17.2222 
16.6666 
15.5555 
16.1111 
17.2222 
16.1111 
16.6666 
18.3333 

19.4444 
16.6666 
16.6666 



RUN NAME 



1 1 14s1_A.PRN 

ST1_3K.PRN 

ST1_30K.PRN 

ST1_50K.PRN 

1121SJB.PRN 

1121SJD.PRN 

1 122S_2A.PRN 

ST2AJ3K.PRN 

ST2A_30K.PRN 

1123SJ3.PRN 

ST37J3K.PRN 

ST37J30K.PRN 

1 123S_7.PRN 

ST73_3K.PRN 

ST73J30K.PRN 

ST7J3K.PRN 

ST7J30K.PRN 

ST10_EDG.PRN 

ST10JDUT.PRN 

1127S_14.PRN 

1128SJ 5.PRN 

1128S_17.PRN 

ST17JN.PRN 

ST17_OUT.PRN 

1128S_18.PRN 

1128SJ9.PRN 

ST19 IN.PRN 

ST19 OUT.PRN 



1 = 1 

PRESSURE RATIO = 0.02968 MACH NUMBER PARAMETER = 0.3741 E-01 

RUN NAME = 111 4s1 A.PRN 



BEGIN NEWTON ITERATION 



ITERATION NUMBER 1 
ITERATION NUMBER 2 
ITERATION NUMBER 3 
ITERATION NUMBER 4 



MACH NO. PARAM. = 0.105475 
MACH NO. PARAM. = 0.104822 
MACH NO. PARAM. = 0.104821 
MACH NO. PARAM. = 0.104821 



ERROR TERM = 
ERROR TERM = 
ERROR TERM = 
ERROR TERM = 



0.6522E-03 

0.1641E-05 

-0.3698E-09 

0.8591E-13 



VREF = 80.23415022742 



I = 28 

PRESSURE RATIO = 0.03004 MACH NUMBER PARAMETER = 0.3787E-01 

RUN NAME = ST19JDUT.PRN 

BEGIN NEWTON ITERATION 



ITERATION NUMBER 
ITERATION NUMBER 
ITERATION NUMBER 



1 MACH NO. PARAM. 

2 MACH NO. PARAM. 

3 MACH NO. PARAM. 



= 0.105777 ERROR 
= 0.105485 ERROR 
= 0.105485 ERROR 



TERM = 0.2921 E-03 
TERM = 0.2872E-06 
TERM = -0.6810E-10 



VREF = 80.51152210108 
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FILE NAME 



EXPERIMENT 

NUMBER 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 



REFERENCE 

VELOCITY 

80.2342 

79.7607 

79.9273 

79.9146 

80.2627 

80.6606 

79.8038 

80.2216 

80.5328 

79.5598 

80.0682 

80.3411 

80.4184 

80.2481 

80.4180 

80.0935 

79.5615 

80.5887 

80.3718 

80.1937 

80.0501 

80.3705 

79.7917 

79.8682 

79.8265 

80.3422 

81.0107 

80.5115 



1114s1_A.PRN 
ST1_3K.PRN 
ST1_30K.PRN 
ST1_50K.PRN 
1121S_1B.PRN 
1121S_1D.PRN 
1122S_2A.PRN 
ST2A_3K.PRN 
ST2A_30K.PRN 
1123S_3.PRN 
ST37_3K.PRN 
ST37_30K.PRN 
1123S_7.PRN 
ST73_3K.PRN 
ST73_30K.PRN 
ST7_3K.PRN 
ST7.30K.PRN 
ST 1 0_EDG.PRN 
ST 1 0_OUT.PRN 
1127S_14.PRN 
1128S_15.PRN 
1128S_17.PRN 
ST17JN.PRN 
ST17_OUT.PRN 
1128S_18.PRN 
1128S_19.PRN 
ST19JN.PRN 
ST19 OUT.PRN 
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APPENDIX C: VELOCITY HISTOGRAMS, TIME-HISTORY PLOTS, POWER- 
SPECTRUM PLOTS, AND AUTO-CORRELATION PLOTS 

1. Station 3, Point 7 



a. 3K Data Points 



Filename* c iNld v50a Ju\st3Nst37_3k . s01 II Processors’ 2 Mode* Coincidence 

10.000 —I 1 i 1 L 



4* 

C 

0 

o 

u 

V 

eu 



7.500 - 



5.000 - 



2.500 - 



0.000 
-238 




903 



-137.901 -36.898 64.104 165.107 

Uelocity (m/s) for Component U 



266.1 



Position (rec/in) (-0.8394, 0.0001, -4.5420) Uelocity at cursor = 0.3 

Uelocity mean = 13.603 Percent at cursor = 2.73 




Position (rec/in) (-0.8394, 0.0001, -4.5420) Uelocity at cursor = 0.4 

Uelocity mean = 16.030 Percent at cursor = 1.81 
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Pos. : (-0.8394, 0.0001, -4.5420) st37_3k.v01 

V Component 



Edit Zoom Mel 



112.972 



-80.912 
















O.00OE+00 



Total Time (ms) 
Data Rate (KHz) 
Mean Velocity 
Standard Dev. 
Vel. at Cursor 
Delta Time (ms) 
4.847E+04 Time at Cursor 



51222 

0.040 

16.03 

32.31 

22.95 

175.0 

25375 



V Component 




0) 

z> 



Data Rate (KHz) 
Mean Velocity 
Standard Dev. 
Vel. at Cursor 



0.040 

13.60 

25.25 

11.10 



Pos. : (-0.8394, 0.0001, -4.5420) st37_3k.v01 I Edit Zoom Hel 

V Component _ < 

Total Time (ms) = 51222 
Data Rate (KHz) = 0.040 
Mean Velocity = 16.03 
Standard Dev. = 32.31 
Vel. at Cursor = -10.0 

. Delta Time (ms) = 5.000 

0 .000E+00 ms 1.385E+03 Time at Cursor = 725.0 



i-n ft— TnJb 


nMli: " F" 


Lr U u u 


nHJ Li u U U _ J 



M 

\ 

£ 



••H 

0 

0 

0) 

D 



89.355 



-62.149 

0.000E+00 ms 1.385E+03 



U Component 



1 ft fin ft flfl- 


n frjift — f 


lr U l/U U u ^ 


vAJIr-’Hk^ f 



Data Rate (KHz) 
Mean Velocity 
Standard Dev. 
Vel. at Cursor 



0.040 

13.60 

25.25 

6.459 
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POWER SPECTRUM POWER SPECTRUM 




FREQUENCV (Hz) 



Graph Component : U 
X of Bins Containing Data = 100 
Scaling Factors (X,¥) = 1.0 , 1.0 




Graph Component : V 
X of Bins Containing Data = 100 
Scaling Factors (X,¥) = 1.0 , 1-0 
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CORRELATION CORRELATION 



i.ODO 



Fi le Tc ' S Id v50a jw\st3Sst37_3k .001 Tit le ’Correlat ion GrapJi 




TIME DELAY (Sec) 



Graph Component I U 

V. of Bins Containing Data = 100 

Scaling Factors (X,Y) = 1.0 , 1.0 




X of Bins Containing Data = 100 
Scaling Factors (X,Y) = 1.0 , 1.0 
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b. 30K Data Points 



Filename^ c :N Id u50a juNst3\st37_30k . sGl if Processors •’ 2 Mode: Coincidence 
10.000 



c 

0 

0 

u 

0 

CL, 



7.500 - 



5.000 - 



2.500 - 



0.000 

-242 




-138.818 -35.633 67.552 170.737 

Velocity (m/s) for Component M 



Position (rec/in) (-0.8393, 0.0000, -4.5420) Velocity at cursor = 0.2 

Velocity mean = 15.972 Percent at cursor = 2.64 



Filename', c :Xldu50a ju\st3\st37_30k .s01 ff Processors: 2 
10.000 -1 1 1 L 



Mode: Coincidence 
L 



P 

0) 

o 

u 

0) 

PU 



7.500 - 



5.000 - 



2.500 - 



0.000 

-310. 




658 -178.881 -47.105 84.672 216.449 

Velocity (m/s) for Component V 



348.2 



Position (rec/in) (—0.8393, 0.0000, —4.5420) Velocity at cursor = 0.1 

Velocity mean = 18.801 Percent at cursor = 1.80 
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Pos. : (-0.8393, 0.0000, -4.5420) st37_30k.v01 

V Component 

7 117.628 



-80.052 







I.U- .iIJ1ji.Ii /hi t K .L. 1. Ahn'.A. 




r rjV ' o iVV r 









0 . OO0E+0O 
93.371 



w 

\ 

£ 



* -61.434 



ms 

M Component 



1.939E+05 









/LA*A I<W Mt ja. 











0 

o 



0 

:> 



0 .000E+00 



ms 



1.939E+05 



| Edit Zoom He 1 

Total Time (ms) = 60012 
Data Rate (KHz) = 0.049 
Mean Velocity = 18.78 
Standard Dev. = 32.94 
Uel . at Cursor = 25.38 
Delta Time (ms) = 700. G 
Time at Cursor = 10150 



Data Rate (KHz) = 0.049 
Mean Velocity = 15.96 
Standard Dev. = 25.80 
Uel. at Cursor = 19.03 



Pos. : (-0.8393, 0.0000, -4.5420) st37_30k . vOl 
7 117.628 



-80.052 



V Component 



n n n r 


rHn il 0 n 


TJYn. fLluWlI U 


nJ n ' JuL 


1_L- j ST 


J — J 1 [_\ U u IQx- 



0.000E+OO ms 1.385E+03 

M Component 



93.371 



w 

\ 

£ 



$ -61.434 






n rrufe 



o 

o 



© 

:> 



"D .000E+00 



ms 



1 .385E+03 



| Edit Zoom Hel 

Total Time (ms) = 60012 
Data Rate (KHz) = 0.049 
Mean Velocity = 18.78 
Standard Dev. = 32.94 
Vel. at Cursor = 73.08 
Delta Time (ms) = 5.00G 
Time at Cursor = 725. G 



Data Rate (KHz) = 0.049 
Mean Velocity = 15.96 
Standard Dev. = 25. 8G 
Vel. at Cursor = 55.25 
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POUER SPECTRUM POWER SPECTRUM 




FREQUENCV (Hz) 



Graph Component : U 

V. of Bins Containing Data = 100 

Scaling Factors (X,V) = 1.0 , 1.0 




FREQUENCV (Hz) 



Graph Component ! V 
X of Bins Containing Data = 100 
Scaling Factors (X,V) = 1.0 , 1.0 
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CORRELATION CORRELATION 



Fi Ie *c *\ldv5Ga ju\st3\st37_30k .UOl Tit le ’Correlation Graph 




X of Bins Containing Data = 100 
Scaling Factors <X,Y) = 1.0 , 1.0 




Graph Component I ^ 

X of Bins Containing Data = 100 
Scaling Factors (X,Y) = 1.0 , 1.0 



63 



2 



Station 7, Point 3 



cl 3K Data Points 




Position (rec/in) (-0.1225, 0.0000, -3.2921) Uelocity at cursor = 19.6 

Uelocity mean = 14.705 Percent at cursor = 1.63 



Filename^ c :Sldv50a JwSst7Sst73_3k .s01 if Processors: 2 Mode: Coincidence 




Position (rec/in) (-0.1225, 0.0000, -3.2921) Uelocity at cursor = 50.0 

Uelocity mean = 39.875 Percent at cursor = 2.57 
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Pos. : (-0.1225, 0.0000, -3.2921) st73_3k.v01 

V Component 




48.256 -I 



U Component 



M 






$ -18.846 



u 

o 

13 



0 . 000E+00 



1 .524E+04 



Edit Zoom Hel 



Total Time (ms) = 37455 
Data Rate (KHz) = 0.082 
Mean Velocity = 39.87 
Standard Dev . = 22. 9S 
Vel. at Cursor = 42.06 
Delta Time (ms) = 55.00 
Time at Cursor = 7975. 



Data Rate (KHz) = 0.002 
Mean Velocity = 14.70 
Standard Dev. = 11.18 
Vel. at Cursor = 16.44 



Pos. : (-0.1225, 0.0000, -3.2921) st73_3k.v01 

V Component 

108.873 



-29 . 123 






0.O0OE+O0 
48.256 



ms 

M Component 



1 .385E+03 



* 

w 

2 -18.846 



'ttt .ifc Tftm 


r/^i , .nnljrmf \ f 




lyi ift/Ms 



u 

o 

V 



0 . 00OE+OO 



ms 



1 . 385E+03 



j Edit Zoom Hel 

Total Time (ms) = 37455 
Data Rate (KHz) = 0.082 
Mean Velocity = 39.87 
Standard Dev. = 22.99 
Vel. at Cursor = 37.21 
Delta Time (ms) = 5.000 
Time at Cursor = 725.0 



Data Rate (KHz) = 0.082 
Mean Velocity = 14.70 
Standard Dev. = 11.18 
Vel. at Cursor = 25.64 
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POWER SPECTRUM POWER SPECTRUM 




FREQUENCY (Hz) 



Graph Component : U 
X of Bins Containing Data = 100 
Scaling Factors (X,V) = 1.0 , 1.0 




FREQUENCV (Hz) 



Graph Component l V 
x of Bins Containing Data = 100 
Scaling Factors (X,Y) = 1.0 f 1.0 
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CORRELATION CORRELATION 




Graph Component : U 
X of Bins Containing Data = 100 
Scaling Factors (X,Y) = 1.0 , 1.0 




Graph Component: V 
X of Bins Containing Data = 100 
Scaling Factors <X,Y) = 1.0 , 1.0 
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b. 



3 OK Data Points 



Filename’ c !\ Idu50a ju\st7Nst73_30k . sOl ft Processors! 2 Mode! Coincidence 
10.000 



7.500 

+> 

P 

0 ) 

2 5.000 

QJ 

o. 



2.500 



0.000 

-103.824 -56.709 -9.595 37.519 84.633 131.7 

Velocity (m/s) for Component M 

Position (rec/in) (-0.1224, 0.0001, -3.2920) Velocity at cursor = 19.1 

Velocity mean = 13.968 Percent at cursor = 1.45 




Filename! c !\ldv50a JuSst7\st73_30k .s01 ft Processors! 2 Mode! Coincidence 



+> 

P 

D 

U 

U 

0 

a, 



10.000 



7.500 - 



5.000 



2.500 - 



0.000 




-208.597 



-110.138 -11.679 86.780 185.240 

Velocity (m/s) for Component \J 



Position (rec/in) (-0.1224, 0.0001, -3.2920) 
Velocity mean = 37.561 



Velocity at cursor 
Percent at cursor = 



283. G 



48.3 

2.40 
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Velocity (m/s) * S’ Velocity (m/s) 



Pos. : (-0.1224, 0.0001, -3.2920) st73_30k.v01 

V Component 

9 111.394 



Edit Zoom Mel 



-36.282 
















































0 .00OE+0O 



ms 

U Component 



Total Time (ms) = 24086 
Data Rate (KHz) = 0.128 
Mean Velocity = 37.55 
Standard Dev. = 24.61 
Vel . at Cursor = 35.01 
Delta Time (ms) = 800.0 
2.216E+05 Time at Cursor = 11600 




Data Rate (KHz) = 0.128 
Mean Velocity = 13.96 
Standard Dev. = 11.77 
Vel. at Cursor = 13.57 



s. : (-0.1224, 0.0001, -3.2920) st73_30k.v01 




0.000E+00 ms 

M Component 

49.299 



1 .385E+03 



-21.368 






j Edit Zoom Hel 

Total Time (ms) = 24086 
Data Rate (KHz) = 0.128 
Mean Velocity = 37.55 
Standard Dev. = 24.61 
Vel. at Cursor = 45.55 
Delta Time (ms) = 5.000 
Time at Cursor = 725.0 



Data Rate (KHz) = 0.128 
Mean Velocity = 13.96 
Standard Dev. = 11.77 
Vel. at Cursor = 19.61 



0 .0OOE+O0 



1 .385E+03 
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POWER SPECTRUM POWER SPECTRUM 




Graph Component : U 

*/ of Bins Containing Data = 100 

Scaling Factors (X,Y) = 1.0 , 1.0 



F i 1 e : c : SI d v50a jwSs t7Ss t73_30k . U0 1 Title: Power Spectrum 



10 u 


2 


-1 




10 A 


~ 


10 -2 


— 


-3 


- 


10 3 


— = 


. -4 


. 


10 * 




-5 


- 


10 D 




IQ" 6 




ID" 7 


’ 




10 



-2 



—i 1 r— 



-i 1 r— 



10 A 10 

FREQUENCY (Hz) 

Graph Component : V 
x of Bins Containing Data = 100 
Scaling Factors (X,Y) = 1.0 , 1.0 



10 * 
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CORRELATION CORRELATION 




X of Bins Containing Data = 100 
Scaling Factors (X,Y) = 1.0 , 1.0 




Graph Component : ^ 

X of Bins Containing Data = 100 
Scaling Factors <X,Y) = 1.0 , 1.0 
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c. 



Ganaim 3K Data Points 



Filename: c:\ldv_bert\st7\0507s71.s03 ft Processors: 2 Mode: Coincidence 



+> 

p 

flj 

o 

u 

0) 

cu 



10. GOG 



7. 500 - 



5. 000 — 



2.500 - 



0.000 

— 107 . 




836 



-58.575 -9.314 39.946 89.207 

Velocity (m/s) for Component U 



138.4 



Position (rec/in) (-0.1224, 0.0000, -3.2921) Velocity at cursor = 20.7 

Velocity mean = 15.316 Percent at cursor = 1.20 




Position (rec/in) (-0.1224, 0.0000, -3.2921) Velocity at cursor = 50.0 

Velocity mean = 38.608 Percent at cursor = 2.41 
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Pos. : (-0.1224, 0.0000, -3.2921) 0507s71 . v03 

V Component 

116.958 



-39.742 



0.000E+00 
52.258 







WflilffMliftiiMflRifll Ifffl 


Jr i H 1T ' "Hn ' » 





ms 

U Component 



j Edit Zoom He! 

Total Time (ms) = 5064. 
Data Rate (KHz) = 0.60G 
Mean Uelocity = 38.60 
Standard Deu. = 26.11 
Uel. at Cursor = 0.113 
Delta Time (ms) = 15.00 
4.155E+03 Time at Cursor = 2175. 



w 

\ 



$ -21 
•H 
0 
o 

p-H 

0) 

D 



633 
0 . O00E+OO 



1 


. -1 . 


i — 




.nw?RfinirciTvmif 

ummmvi'ifmmrf 


|* 1 H 1 | ' Ij • n * ** 







Data Rate (KHz) = 0.607 
Mean Uelocity = 15.31 
Standard Dev. = 12.31 
Uel. at Cursor = -2.67 
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Edit Zoom Hel 



Total Time (ms) = 5064. 
Data Rate (KHz) = 0.606 
Mean Uelocity = 38.60 
Standard Dev. = 26.11 
Uel. at Cursor = 53.25 
Delta Time (ms) = 5.000 
Time at Cursor = 725.0 



Data Rate (KHz) = 0.607 
Mean Uelocity = 15.31 
Standard Dev. = 12.31 
Uel. at Cursor = 20.37 
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POWER SPECTRUM POWER SPECTRUM 



FI le IctSldu hert\st7\0SD7s71 T I t 1 I Pnupr' Rnpr+rtim 




FREQUENCY (Hz) 



Graph Component : U 
X of Bins Containing Data = 100 
Scaling Factors (X,Y) = 1.0 , 1.0 




FREQUENCY (Hz) 



Graph Component l V 

V. of Bins Containing Data = 100 

Scaling Factors (X,Y) = 1.0 , 1.0 
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CORRELATION CORRELATION 




X of Bins Containing Data = 100 
Scaling Factors (X,Y) = 1.0 , 1.0 



File Tc *N Id v bert\st7\0507s71 .U03 Tit le :Correlat ion Graph 




X of Bins Containing Data = 100 
Scaling Factors (X,Y) = 1.0 , 1.0 
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3. 



Station 7, Point 20 



a. 3K Data Points 




Position (rec/in) (0.3654, 0.0000, -3.2920) Velocity at cursor = 30.9 

Velocity mean = 26.004 Percent at cursor = 3.22 




Position (rec/in) (0.3654, 0.0000, -3.2920) Velocity at cursor = 68.9 

Velocity mean = 59.130 Percent at cursor = 1.73 
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Pos. : (0.3654, 0.0000, -3.2920) st7_3k.u01 

V Component 

126.494 



-8.275 
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4 .432E+04 



Edit Zoom He 1 



Total Time (ms) = 44037 
Data Rate (KHz) = 0.070 
Mean Velocity = 59.10 
Standard Deu. = 22.46 
Vel. at Cursor = 68.80 
Delta Time (ms) = 160.0 
Time at Cursor = 23200 



Data Rate (KHz) = 0.070 
Mean Velocity = 26.80 
Standard Deu. = 9.219 
Vel. at Cursor = 29.98 



Pos. : (0.3654, 0.0000, -3.2920) st7_3k.u01 

\t Component 



126.494 



-8.275 
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U Component 



1.385E+03 



| Edit Zoom Hel 

Total Time (ms) = 44037 
Data Rate (KHz) = 0.070 
Mean Velocity = 59.10 
Standard Deu. = 22.46 
Vel. at Cursor = 74.07 
Delta Time (ms) = 5.000 
Time at Cursor = 725.0 
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Data Rate (KHz) 
Mean Velocity 
Standard Deu . 
Vel. at Cursor 



0.070 

26.88 

9.219 

28.07 
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POUER SPECTRUM POWER SPECTRUM 




FREQUENCY (Hz) 



Graph Component : U 
X of Bins Containing Data = IDG 
Scaling Factors (X,Y) = 1.0 , 1.0 



File:c:\ldu50a iw\st7Sst7 3k.V01 Title :Power Spectrum 




Graph Component : V 
X of Bins Containing Data - 100 
Scaling Factors (X,Y) = l.G , l.G 
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CORRELATION CORRELATION 




V. of Bins Containing Data = 100 
Scaling Factors (X,Y) = 1.0 , 1.0 



1.000 



0.659 



CT.317 



-0.024 



-0.366 



Fi Ie Tc :\ldv50a ju\st7\st7_3k .U01 Ti t Ie :Corre lat Ion Graph 

I t I I I I 1 I I I I I 1 I — 1 1 1 1 L 




Time Delay (sec) 



Graph Component : V 

y. of Bins Containing Data = 100 

Scaling Factors (X,Y) = 1.0 , 1.0 
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b. 



3 OK Data Points 




Position (rec/in) (0.3653, -0.0001, —3.2920) Velocity at cursor = 30. G 

Velocity mean = 25.181 Percent at cursor = 6.04 




Position (rec/in) (0.3653, -0.0001, -3.2920) Velocity at cursor = 64.5 

Velocity mean = 53.312 Percent at cursor = 1.83 
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Pos. 

2 129.733 



(0.3653, -0.0001, -3.2920) st7__30k . uOl 
V Component 



-23.117 
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-8 . 149 



ms 

U Component 



I Edit Zoom Hel 

Total Time (ms) = 43243 
Data Rate (KHz) = 0.071 
Mean Velocity = 53. 3G 
Standard Dev . = 25.47 

Vel. at Cursor = 56.62 
Delta Time (ms) = 1000. 
2.770E+05 Time at Cursor = 14500 




















0.000E+00 



Data Rate (KHz) = 0.071 
Mean Velocity = 25.18 
Standard Dev. = 11.11 
Vel. at Cursor = 25.27 



2 .770E+05 



Pos. : (0.3653, -0.0001, -3.2920) st7_30k . vOl 

V Component 



Edit Zoom Hel 



2 129.733 



-23.117 
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T3.000E+00 



Total Time (ms) = 43243 
Data Rate (KHz) = 0.071 
Mean Velocity = 53.30 
Standard Dev. = 25.47 
Vel. at Cursor = 60.62 
Delta Time (ms) = 5.000 
1.385E+03 Time at Cursor = 725.0 




Data Rate (KHz) = 0.071 
Mean Velocity = 25.18 
Standard Dev. = 11.11 
Vel. at Cursor = 31.40 
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POWER SPECTRUM 




Graph Component : U 
X of Bins Containing Data = 100 
Scaling Factors (X,Y) = 1.0 , 1.0 




FREQUENCY (Hz) 



Graph Component * V 
X of Bins Containing Data = 100 
Scaling Factors (X,Y) = 1.0 , 1.0 
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CORRELATION 




y. o£ Bins Containing Data = 100 
Scaling Factors (X,Y) = 1.0 , 1.0 




y. of Bins Containing Data = 100 
Scaling Factors (X,Y) = 1.0 , 1.0 
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c. 



Ganaim 3K Data Points 




Position (rec/in) (0.3654, 0.0000, -3.2921) Velocity at cursor = 31.4 

Velocity mean = 26.554 Percent at cursor = 4.82 




Position (rec/ln) (0.3654, 0.0000, -3.2921) Velocity at cursor = 68.7 

Velocity mean = 57.071 Percent at cursor = 2.60 
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Pos. : (0.3654, 0.0000, -3.2921) 0507s71 . u20 

V Component 

136.833 



-22.691 
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U Component 
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Edit Zoom He 1 



Total Time (ms) = 5012. 
Data Rate (KHz) = 0.613 
Mean Uelocity = 57.07 
Standard Dev. = 26.50 
Uel. at Cursor = 39.96 
Delta Time (ms) = 20.00 
Time at Cursor = 2900. 



Data Rate (KHz) = 0.613 
Mean Uelocity = 26.55 
Standard Deu. = 11.08 
Uel. at Cursor = 24.18 
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0 . 00OE+00 



ms 



Edit Zoom Hel 



Total Time (ms) = 5012. 
Data Rate (KHz) = 0.613 
Mean Uelocity = 57.07 
Standard Deu. = 26.58 
Uel. at Cursor = 77.97 
Delta Time (ms) = 5.000 
Time at Cursor = 725.0 



Data Rate (KHz) = 0.613 
Mean Uelocity = 26.55 
Standard Deu. = 11.08 
Uel. at Cursor = 29.97 



1 .385E+03 
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POWER SPECTRUM POWER SPECTRUM 




Graph Component I U 

y, of Bins Containing Data = 100 

Scaling Factors (X,Y) = 1.0 , 1.0 




FREQUENCY (Hz) 



Graph Component: V 

y, of Bins Containing Data = 100 

Scaling Factors (X,Y) = 1.0 , 1.0 
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CORRELATION CORRELATION 



1 .OGO 



G .700 



0.401 



0.101 



-0 . 198 



F i Ie :c T\ldv_ber t\st7\0507s71 .U20 T it le iCorrelat ion Graph 




Time Delay (sec) 



Graph Component : U 
X of Bins Containing Data = 100 
Scaling Factors (X,V) = 1.0 , 1.0 



FilercAIdv bert\st7S0507s71 .U20 Tit le :Correlat ion Graph 




X of Bins Containing Data = 100 
Scaling Factors (X,V) = 1.0 , 1.0 
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APPENDIX D: REDUCED DATA FROM STATIONS 1, 3, 7, 19 



Pitchwise Survey at Station 1 



X (in) 


Y (in) 


U/Vref 


V/Vref 


UVMag/Vref 


U - Turbulence 


V - Turbulence 


Reynolds 

Stress 


Correlation 

Coefficient 


1.00 


-6.29 


0.7560 


0.6273 


0.9824 


4.5167 


4.2288 


-0.5645 


-0.0459 


0.75 


-6.29 


0.7527 


0.6348 


0.9847 


4.2846 


4.1545 


-0.1621 


-0.0141 


0.50 


-6.29 


0.7549 


0.6361 


0.9872 


4.1377 


4.1513 


-0.8599 


-0.0778 


0.25 


-6.29 


0.7532 


0.6393 


0.9879 


3.9633 


4.2387 


-0.0287 


-0.0027 


0.00 


-6.29 


0.7551 


0.6417 


0.9910 


4.2522 


4.0219 


-0.0638 


-0.0058 


-0.25 


-6.29 


0.7573 


0.6399 


0.9914 


4.1720 


4.0454 


-0.0756 


-0.0070 


-0.50 


-6.29 


0.7646 


0.6367 


0.9950 


4.3993 


3.9576 


0.4649 


0.0415 


-0.75 


-6.29 


0.7671 


0.6362 


0.9966 


4.3408 


3.9461 


0.1865 


0.0169 


-1.00 


-6.29 


0.7655 


0.6322 


0.9929 


4.5415 


3.8821 


0.4439 


0.0391 


-1.25 


-6.29 


0.7685 


0.6303 


0.9940 


4.0621 


3.6989 


0.3160 


0.0327 


-1.50 


-6.29 


0.7622 


0.6312 


0.9896 


4.9631 


3.7591 


0.0913 


0.0076 


-1.75 


-6.29 


0.7527 


0.6333 


0.9837 


5.6056 


3.5044 


0.2665 


0.0211 


-2.00 


-6.29 


0.7505 


0.6371 


0.9844 


4.0228 


3.4712 


-0.0511' 


-0.0057 


-2.25 


-6.29 


0.7509 


0.6368 


0.9846 


3.6099 


3.5205 


-0.0113 


-0.0014 


-2.50 


-6.29 


0.7519 


0.6417 


0.9885 


3.5904 


3.3798 


0.1693 


0.0217 


-2.75 


-6.29 


0.7512 


0.6462 


0.9909 


3.5229 


3.4416 


-0.1665 


-0.0213 


-3.00 


-6.29 


0.7553 


0.6495 


0.9961 


3.5721 


3.2735 


0.1392 


0.0185 


-3.25 


-6.29 


0.7576 


0.6500 


0.9982 


3.9669 


3.3235 


-0.2241 


-0.0264 


-3.50 


-6.29 


0.7767 


0.6438 


1.0088 


3.5708 


3.4123 


-0.3948 


-0.0503 


-3.75 


-6.29 


0.7787 


0.6344 


1.0044 


3.5851 


3.5118 


-0.3790 


-0.0468 


-4.00 


-6.29 


0.7753 


0.6322 


1.0004 


3.4988 


3.3147 


-0.1981 


-0.0265 


-4.25 


-6.29 


0.7714 


0.6303 


0.9962 


3.4435 


3.3032 


-0.1314 


-0.0179 


-4.50 


-6.29 


0.7634 


0.6295 


0.9894 


3.7142 


3.1308 


0.0422 


0.0056 


-4.75 


-6.29 


0.7596 


0.6341 


0.9895 


3.4552 


2.9855 


0.2140 


0.0322 


-5.00 


-6.29 


0.7592 


0.6383 


0.9918 


3.3208 


2.9673 


-0.2577 


-0.0406 


-5.25 


-6.29 


0.7596 


0.6459 


0.9971 


2.9610 


2.8494 


-0.2546 


-0.0469 


-5.50 


-6.29 


0.7652 


0.6516 


1.0050 


3.2893 


2.7745 


0.0387 


0.0066 


-5.75 


-6.29 


0.7714 


0.6537 


1.0111 


2.9942 


3.0307 


-0.3735 


-0.0639 


-6.00 


-6.29 


0.7748 


0.6563 


1.0154 


3.2675 


2.9769 


-0.1074 


-0.0171 


-6.25 


-6.29 


0.7782 


0.6528 


1.0157 


2.9052 


2.6739 


-0.1011 


-0.0202 


-6.50 


-6.29 


0.7787 


0.6490 


1.0137 


3.4300 


2.7946 


-0.0564 


-0.0091 


-6.75 


-6.29 


0.7773 


0.6426 


1.0085 


2.7609 


2.5602 


-0.0539 


-0.0119 


-7.00 


-6.29 


0.7767 


0.6388 


1.0057 


3.0063 


2.5740 


0.4140 


0.0831 


-7.25 


-6.29 


0.7741 


0.6370 


1.0024 


3.1116 


2.4751 


0.2009 


0.0405 


-7.50 


-6.29 


0.7727 


0.6406 


1.0038 


2.7515 


2.3997 


0.2478 


0.0583 


-7.75 


-6.29 


0.7686 


0.6405 


1.0005 


2.5751 


2.2243 


0.1246 


0.0338 


-8.00 


-6.29 


0.7665 


0.6444 


1.0014 


2.8413 


2.2209 


0.1098 


0.0270 


-8.25 


-6.29 


0.7686 


0.6469 


1.0046 


3.1292 


2.2872 


0.0908 


0.0197 


-8.50 


-6.29 


0.7694 


0.6525 


1.0089 


3.0650 


2.0272 


0.0634 


0.0158 


-8.75 


-6.29 


0.7739 


0.6534 


1.0129 


2.9483 


2.1966 


0.0291 


0.0070 


-9.00 


-6.29 


0.7783 


0.6542 


1.0167 


2.4833 


2.1338 


0.0240 


0.0070 
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Pitchwise Survey at Station 3 





Total Flow Average 






Negative Average 


Positive Average 


X (in) 


Y<«) 


U/Vref 


V/Vref 


Utot/Vref 


U-Turb 


V-Turb 


U/Vref 


V/Vref 


U/Vref 


V/Vref 


-0.9160 


-4.5421 


-0.0508 


-0.0820 


0.0964 


-6.8660 


-7.5121 


-0.0745 


-0.1005 


0.0976 


0.1242 


-0.9061 


-4.5420 


-0.0566 


-0.0813 


0.0991 


-7.0072 


-9.0429 


-0.0780 


-0.1061 


0.0952 


0.1348 


-0.8951 


-4.5420 


-0.0566 


-0.0807 


0.0986 


-7.3448 


-10.6218 


-0.0814 


-0.1091 


0.1169 


0.2125 


-0.8831 


-4.5420 


-0.0524 


-0.0713 


0.0885 


-10.4817 


-13.7893 


-0.0855 


-0.1150 


0.1897 


0.2506 


-0.8699 


-4.5420 


-0.0105 


-0.0184 


0.0212 


-18.3077 


-24.1420 


-0.0848 


-0.1154 


0.3524 


0.4773 


-0.8553 


-4.5420 


0.0692 


0.0767 


0.1033 


26.3741 


33.9905 


-0.0825 


-0.1146 


0.4479 


0.5715 


-0.8393 


-4.5420 


0.2345 


0.2896 


0.3726 


33.1398 


43.3500 


-0.0812 


-0.1178 


0.5142 


0.6610 


-0.8217 


-4.5420 


0.3414 


0.4153 


0.5377 


33.9309 


43.6112 


-0.0799 


-0.1182 


0.5525 


0.6905 


-0.8023 


-4.5420 


0.4373 


0.5191 


0.6787 


33.0556 


41.3025 


-0.0764 


-0.1118 


0.5858 


0.7099 


-0.7811 


-4.5420 


0.5376 


0.6040 


0.8086 


29.6847 


36.5104 


-0.0753 


-0.1088 


0.6291 


0.7162 


-0.7575 


-4.5420 


0.5532 


0.5886 


0.8078 


28.0982 


34.1066 


-0.0704 


-0.1176 


0.6245 


0.6757 


-0.7317 


-4.5420 


0.6236 


0.6658 


0.9122 


23.7893 


30.3195 


-0.0718 


-0.0978 


0.6634 


0.7168 


-0.7035 


-4.5420 


0.6598 


0.6660 


0.9375 


20.8838 


26.3703 


-0.0571 


-0.0992 


0.6841 


0.6953 


-0.6723 


-4.5420 


0.6963 


0.6885 


0.9792 


17.4786 


23.8027 


-0.0673 


-0.1251 


0.7080 


0.7041 


-0.6377 


-4.5420 


0.7192 


0.7057 


1.0076 


14.4658 


22.1652 


-0.0455 


-0.1565 


0.7255 


0.7169 


-0.6001 


-4.5420 


0.7387 


0.7207 


1.0320 


12.5289 


20.1208 


0.0000 


0.0000 


0.7422 


0.7275 


-0.5585 


-4.5420 


0.7553 


0.7316 


1.0515 


10.5465 


18.8519 


0.0000 


0.0000 


0.7565 


0.7352 


-0.5127 


-4.5420 


0.7781 


0.6664 


1.0245 


10.0197 


15.9950 


0.0000 


0.0000 


0.7794 


0.6687 


-0.4626 


-4.5420 


0.7642 


0.7640 


1.0806 


8.0755 


15.7110 


0.0000 


0.0000 


0.7642 


0.7643 


-0.4073 


-4.5420 


0.7645 


0.7642 


1.0809 


8.4760 


15.8142 


0.0000 


0.0000 


0.7645 


0.7647 


-0.3465 


-4.5420 


0.7524 


0.7933 


1.0933 


8.0195 


15.1845 


0.0000 


0.0000 


0.7524 


0.7933 


-0.2796 


-4.5420 


0.7442 


0.7777 


1.0764 


8.2561 


14.5299 


0.0000 


0.0000 


0.7442 


0.7777 


-0.2060 


-4.5420 


0.7460 


0.7736 


1.0746 


8.6669 


13.4328 


0.0000 


0.0000 


0.7460 


0.7736 


-0.1251 


-4.5420 


0.7330 


0.7705 


1.0635 


9.0455 


13.2214 


0.0000 


0.0000 


0.7330 


0.7705 


-0.0361 


-4.5420 


0.7322 


0.7560 


1.0525 


9.3525 


12.0596 


0.0000 


0.0000 


0.7322 


0.7560 


0.0619 


-4.5420 


0.7208 


0.7548 


1.0437 


9.5275 


11.5091 


0.0000 


0.0000 


0.7208 


0.7548 


0.1696 


-4.5420 


0.7186 


0.7403 


1.0317 


9.6073 


10.9092 


0.0000 


0.0000 


0.7186 


0.7406 


0.2882 


-4.5420 


0.6990 


0.7356 


1.0148 


9.5020 


10.1140 


0.0000 


0.0000 


0.6990 


0.7356 


0.4185 


-4.5420 


0.6903 


0.7245 


1.0007 


9.2803 


9.0998 


0.0000 


0.0000 


0.6903 


0.7245 


0.5620 


-4.5420 


0.6706 


0.7162 


0.9811 


9.0082 


8.4316 


0.0000 


0.0000 


0.6706 


0.7162 


0.7196 


-4.5420 


0.6449 


0.6946 


0.9479 


8.7474 


8.0618 


0.0000 


0.0000 


0.6449 


0.6946 


0.8932 


-4.5420 


0.6282 


0.6704 


0.9187 


7.8477 


7.6788 


0.0000 


0.0000 


0.6282 


0.6704 
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Pitchwise Survey at Station 7 





Total Flow Average 






Average of Negative 


Average of Positive 


X(m) 


v (in) 


U/vref 


V/Vref 


Utot/Vxef 


U-Turb 


V-Turb 


U/Vref 


V/Vref 


U/Vref 


V/Vref 


-0.1435 


-3.2920 


0.1749 


0.4748 


0.5060 


12.2241 


26.7540 


-0.0643 


-0.1126 


0.2064 


0.5479 


-0.1335 


-3.2920 


0.1809 


0.4820 


0.5148 


12.9712 


26.9304 


-0.0657 


-0.1081 


0.2146 


0.5524 


-0.1225 


-3.2920 


0.1780 


0.4843 


0.5159 


13.8963 


28.4924 


-0.0656 


-0.1100 


0.2211 


0.5748 


-0.1105 


-3.2920 


0.1861 


0.4919 


0.5259 


14.7300 


30.1213 


-0.0748 


-0.1079 


0.2333 


0.5898 


-0.0973 


-3.2920 


0.1889 


0.4992 


0.5338 


15.0800 


30.5573 


-0.0691 


-0.1077 


0.2385 


0.6005 


-0.0827 


-3.2920 


0.1857 


0.4874 


0.5216 


15.8231 


31.1735 


-0.0719 


-0.1014 


0.2387 


0.5935 


-0.0667 


-3.2920 


0.1938 


0.5026 


0.5387 


15.9031 


32.6549 


-0.0739 


-0.1024 


0.2487 


0.6161 


-0.0481 


-3.2920 


0.2026 


0.5130 


0.5515 


16.1680 


33.0578 


-0.0707 


-0.1090 


0.2588 


0.6283 


-0.0297 


-3.2920 


0.2116 


0.5290 


0.5698 


16.4695 


33.2063 


-0.0736 


-0.0988 


0.2664 


0.6396 


-0.0084 


-3.2920 


0.2166 


0.5227 


0.5658 


16.8288 


35.1634 


-0.0681 


-0.0991 


0.2776 


0.6503 


0.0151 


-3.2920 


0.2136 


0.5201 


0.5622 


16.7788 


36.0857 


-0.0675 


-0.0983 


0.2728 


0.6490 


0.0409 


-3.2920 


0.2162 


0.5202 


0.5634 


16.9932 


36.5286 


-0.0664 


-0.0981 


0.2778 


0.6529 


0.0692 


-3.2920 


0.2401 


0.5600 


0.6093 


17.0275 


37.3337 


-0.0654 


-0.0937 


0.3014 


0.6985 


0.1004 


-3.2920 


0.2453 


0.5712 


0.6216 


16.7892 


37.0429 


-0.0686 


-0.0993 


0.2974 


0.6905 


0.1348 


-3.2920 


0.2533 


0.5765 


0.6297 


16.9775 


38.0062 


-0.0604 


-0.0986 


0.3091 


0.7065 


0.1725 


-3.2920 


0.2580 


0.5909 


0.6448 


16.8549 


37.0968 


-0.0574 


-0.0974 


0.3106 


0.7097 


0.2141 


-3.2920 


0.2733 


0.6101 


0.6685 


16.1963 


36.6200 


-0.0642 


-0.0975 


0.3161 


0.7126 


0.2599 


-3.2920 


0.2849 


0.6250 


0.6869 


15.8988 


36.1124 


-0.0660 


-0.0930 


0.3271 


0.7290 


0.3100 


-3.2920 


0.2988 


0.6631 


0.7273 


15.0613 


33.7251 


-0.0568 


-0.0940 


0.3324 


0.7409 


0.3654 


-3.2920 


0.3156 


0.7053 


0.7726 


14.0393 


31.8136 


-0.0652 


-0.0853 


0.3420 


0.7700 


0.4259 


-3.2920 


0.3282 


0.7261 


0.7969 


12.8077 


29.6575 


-0.0678 


-0.0934 


0.3481 


0.7733 


0.4930 


-3.2920 


0.3323 


0.7293 


0.8014 


12.6994 


27.7864 


-0.0715 


-0.1135 


0.3496 


0.7683 


0.5666 


-3.2920 


0.3469 


0.7369 


0.8145 


11.3452 


25.8805 


-0.0583 


-0.0838 


0.3550 


0.7616 


0.6474 


-3.2920 


0.3549 


0.7283 


0.8102 


11.2177 


24.6712 


-0.0594 


-0.0810 


0.3644 


0.7534 


0.7365 


-3.2920 


0.3667 


0.7457 


0.8310 


10.4527 


22.1508 


-0.0553 


-0.0619 


0.3716 


0.7615 


0.8345 


-3.2920 


0.3819 


0.7624 


0.8527 


9.8618 


18.9547 


0.0000 


0.0000 


0.3835 


0.7707 


0.9422 


-3.2920 


0.3868 


0.7551 


0.8484 


10.0925 


17.5833 


0.0000 


-0.0553 


0.3885 


0.7614 


1.0608 


-3.2920 


0.3912 


0.7689 


0.8627 


9.6454 


15.1138 


0.0000 


0.0000 


0.3914 


0.7730 


1.1911 


-3.2920 


0.3919 


0.7749 


0.8683 


9.1682 


12.5425 


0.0000 


0.0000 


0.3921 


0.7761 


1.3345 


-3,2920 


0.3849 


0.7814 


0.8710 


8.8570 


10.8529 


0.0000 


0.0000 


0.3849 


0.7814 


1.4922 


-3.2920 


0.3739 


0.7863 


0.8706 


7.9523 


8.9572 


0.0000 


0.0000 


0.3739 


0.7863 


1.6657 


-3.2920 


0.3500 


0.7726 


0.8482 


6.8941 


7.5859 


0.0000 


0.0000 


0.3500 


0.7729 
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Total Flow Average 






Average of Negative 


Average of Positive 


X (inches) 


Y (inches) 


U/Vref 


V/Vref 


UVtot/ 

Vref 


U-Turb 


V-Turb 


U/Vref 


V/Vref 


U/Vref 


V/Vref 


3.9999 


1.0620 


0.0253 


0.6604 


0.6609 


7.7800 


24.6470 


-0.0637 


-0.0788 


0.0565 


0.7056 


3.8750 


1.0620 


0.0237 


0.6190 


0.6194 


8.1843 


25.8908 


-0.0709 


-0.0775 


0.0604 


0.6619 


3.7499 


1.0620 


0.0195 


0.5861 


0.5864 


8.5960 


25.6453 


-0.0669 


-0.0846 


0.0684 


0.6254 


3.6250 


1.0620 


0.0197 


0.5217 


0.5220 


9.1367 


24.5519 


-0.0690 


-0.0831 


0.0751 


0.5499 


3.5001 


1.0620 


0.0198 


0.4731 


0.4735 


9.5804 


22.3559 


-0.0723 


-0.0780 


0.0751 


0.4960 


3.3751 


1.0620 


0.0168 


0.4034 


0.4038 


10.0333 


19.5244 


-0.0739 


-0.0786 


0.0777 


0.4141 


3.2499 


1.0620 


0.0312 


0.3255 


0.3270 


12.0674 


17.0258 


-0.0693 


-0.0608 


0.1015 


0.3312 


3.1251 


1.0620 


0.0546 


0.3268 


0.3313 


13.4243 


15.0200 


-0.0629 


-0.0498 


0.1278 


0.3488 


3.0000 


1.0620 


0.0888 


0.5394 


0.5466 


10.7297 


15.4455 


0.0000 


0.0000 


0.1275 


0.5590 


2.8750 


1.0620 


0.0594 


0.7741 


0.7764 


6.3753 


11.0721 


0.0000 


0.0000 


0.0748 


0.7823 


2.7498 


1.0620 


0.0458 


0.8239 


0.8251 


5.0207 


5.1598 


0.0000 


0.0000 


0.0567 


0.8254 


2.6249 


1.0620 


0.0428 


0.8278 


0.8289 


4.3926 


3.6575 


0.0000 


0.0000 


0.0521 


0.8267 


2.5000 


1.0620 


0.0412 


0.8272 


0.8282 


3.9808 


3.3833 


0.0000 


0.0000 


0.0488 


0.8268 


2.375 1 


1.0620 


0.0394 


0.8275 


0.8284 


3.5717 


3.0630 


0.0000 


0.0000 


0.0459 


0.8272 


2.2500 


1.0620 


0.0383 


0.8259 


0.8268 


3.7252 


3.5124 


0.0000 


0.0000 


0.0455 


0.8261 


2.1251 


1.0620 


0.0403 


0.8262 


0.8272 


4.5719 


5.5106 


0.0000 


0.0000 


0.0498 


0.8268 


2.0001 


1.0620 


0.0399 


0.8256 


0.8266 


4.3079 


5.3561 


-0.0098 


-0.1220 


0.0484 


0.8273 


1.8750 


1.0620 


0.0387 


0.8217 


0.8227 


4.5239 


6.3894 


0.0000 


0.0000 


0.0485 


0.8232 


1.7500 


1.0620 


0.0378 


0.8125 


0.8134 


5.1930 


8.7774 


0.0000 


0.0000 


0.0486 


0.8173 


1.6250 


1.0619 


0.0368 


0.8066 


0.8075 


5.1928 


10.1244 


0.0000 


0.0000 


0.0474 


0.8132 


1.5000 


1.0620 


0.0500 


0.7315 


0.7332 


8.6155 


19.3049 


-0.0967 


41.0481 


0.0707 


0.7458 


1.3750 


1.0620 


0.0507 


0.7019 


0.7038 


8.7904 


22.0152 


-0.0643 


41.0590 


0.0711 


0.7229 


1.2500 


1.0620 


0.0410 


0.7076 


0.7087 


8.1040 


22.2274 


-0.0594 


-0.0786 


0.0626 


0.7366 


1.1250 


1.0621 


0.0373 


0.6794 


0.6804 


7.8614 


24.8061 


-0.0612 


-0.0784 


0.0596 


0.7198 


1.0000 


1.0620 


0.0393 


0.6499 


0.6511 


8.0620 


26.0936 


-0.0509 


-0.0969 


0.0640 


0.6913 


0.8750 


1.0620 


0.0328 


0.6525 


0.6533 


7.8186 


24.9827 


-0.0613 


-0.0801 


0.0620 


0.6917 


0.7499 


1.0620 


0.0304 


0.5986 


0.5994 


8.5539 


25.4173 


-0.0544 


-0.0882 


0.0691 


0.6326 


0.6250 


1.0620 


0.0305 


0.5419 


0.5428 


8.3361 


23.6395 


-0.0700 


41.0797 


0.0724 


0.5647 


0.5000 


1.0620 


0.0267 


0.4712 


0.4720 


8.5616 


22.6358 


-0.0644 


-0.0942 


0.0747 


0.4977 


0.3750 


1.0620 


0.0321 


0.3932 


0.3945 


10.3141 


19.8753 


-0.0570 


-0.0763 


0.0825 


0.4085 


0.2500 


1.0620 


0.0439 


0.3235 


0.3265 


12.6135 


16.9451 


-0.0607 


-0.0652 


0.1039 


0.3344 


0.1250 


1.0620 


0.0717 


0.3254 


0.3332 


13.0058 


14.8209 


-0.0643 


-0.0313 


0.1275 


0.3390 


-0.0001 


1.0620 


0.0900 


0.5310 


0.5386 


10.3600 


16.2890 


0.0000 


0.0000 


0.1242 


0.5487 


-0.1250 . 


. 1.0620 


0.0575 


0.7944 


0.7965 


6.4601 


10.7090 


0.0000 


0.0000 


0.0731 


0.8025 


-0.2501 


1.0620 


0.0476 


0.8425 


0.8438 


5.0608 


4.6542 


0.0000 


0.0000 


0.0582 


0.8430 


-0.3751 


1.0620 


0.0417 


0.8439 


0.8449 


4.5586 


3.8154 


0.0000 


0.0000 


0.0520 


0.8431 


-0.4999 


1.0620 


0.0388 


0.8446 


0.8454 


4.2244 


3.5776 


0.0000 


0.0000 


0.0488 


0.8437 


-0.6250 


1.0620 


0.0384 


0.8438 


0.8447 


3.8185 


3.2821 


0.0000 


0.0000 


0.0467 


0.8436 


-0.7500 


1.0620 


0.0365 


0.8424 


0.8432 


3.7316 


3.2875 


0.0000 


0.0000 


0.0449 


0.8419 


-0.8750 


1.0620 


0.0362 


0.8393 


0.8400 


3.6055 


3.7414 


0.0000 


0.0000 


0.0437 


0.8395 


- 1.0000 


1.0620 


0.0340 


0.8357 


0.8364 


3.5738 


4.0021 


0.0000 


0.0000 


0.0418 


0.8367 


-1.1251 


1.0620 


0.0325 


0.8338 


0.8345 


3.6817 


4.7620 


0.0000 


0.0000 


0.0417 


0.8361 


-1.2500 


1.0620 


0.0343 


0.8329 


0.8336 


3.7131 


5.2161 


0.0000 


0.0000 


0.0412 


0.8347 


-1.3750 


1.0620 


0.0323 


0.8272 


0.8279 


3.9602 


6.7284 


0.0000 


0.0000 


0.0398 


0.8296 


-1.5000 


1.0620 


0.0305 


0.8207 


0.8213 


4.4862 


8.7431 


-0.1331 


-0.0914 


0.0401 


0.8274 


-1.6250 


1.0621 


0.0331 


0.8017 


0.8024 


5.1570 


13.3611 


-0.0295 


-0.1048 


0.0441 


0.8119 


-1.7500 


1.0620 


0.0429 


0.7382 


0.7394 


7.8387 


20.1948 


-0.0238 


-0.0818 


0.0646 


0.7614 


-1.8750 


1.0620 


0.0379 


0.7059 


0.7069 


7.8834 


22.3243 


-0.0529 


-0.0771 


0.0637 


0.7381 


-2.0000 


1.0620 


0.0318 


0.6920 


0.6927 


7.6868 


23.7145 


-0.0669 


-0.0788 


0.0598 


0.7296 
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APPENDIX E: CFD GRID GENERATION 



The FORTRAN grid generation program employed was the modified version of 
GRAPE (GRids About Aerofoils using Poisson’s Equation) by R. V. Chima, [Ref. 11]. 
Prior to this study, the grid’s leading edge was blunt or flat-faced, which led to 
inaccuracies in the computed solutions. It was determined that GRAPE properly 
smoothed the points around the leading edge, and that it required definition of the leading 
edge with more user-supplied coordinates. Careful attention was given to properly 
blending the top and bottom grid surfaces into the 0.045 inch radius circular arc that 
defines the leading edge. This was done by ensuring that first derivative continuity was 
within a few percent difference. The final coordinates for the CD blade definition are 
contained below in namelist &GRID3. 



&GRID1 

JMAX=240,KMAX=49,NTETYP=3,NAIRF=5,JAIRF=163,NIBDST=7, 

DSI=0.0002, JTEBOT =60, JTETOP= 181 ,XLE=0.0,N OBSHP=7,XTE= 1 . 0, 
ALAMF=0.0,ALAMR=0.0,XLEFT=-.3,XRIGHT=3., 

RCORN=. 10,NOUT=4,NORDA=5,3,MAXrrA=100, 100, 

&END 

&GRID2 

AAAI=.3,BBBI=.3,DSOBI=0.0 l,ROTANG=-14 .4, 
XTFRAC=1.5,PITCH=.5988,DSRA=.4920,DSLE=.0008,DSTE=.001, 
NLE=28,NTE=18,WAKEP=1 ,OMEGR=1.0,OMEGS=1 0,OMEGP=1 0,OMEGQ=1 .0, 
&END 

&GRID3 AIRFX= 

0.9996, 0.9971, 0.9929, 

0.9601, 0.9453, 0.9305, 

0.8419, 0.8271, 0.8124, 

0.7238, 0.7090, 0.6942, 

0.6056, 0.5908, 0.5760, 

0.4874, 0.4726, 0.4579, 

0.3693, 0.3545, 0.3397, 

0.2511, 0.2364, 0.2216, 

0.1329, 0.1181, 0.1034, 

0.0244, 0.0114, 0.0097, 

0.0033, 0.0021, 0.0011, 

0.0010, 0.0019, 0.0030, 

0.0591, 0.0738, 0.0886, 

0.1772, 0.1920, 0.2068, 

0.2954, 0.3102, 0.3250, 

0.4136, 0.4283, 0.4431, 

0.5317, 0.5465, 0.5612, 



0.9876, 0.9830, 0.9803, 0.9776, 0.9749, 
0.9157, 0.9010, 0.8862, 0.8714, 0.8567, 
0.7976, 0.7828, 0.7681, 0.7533, 0.7385, 
0.6795, 0.6647, 0.6499, 0.6352, 0.6204, 
0.5612, 0.5465, 0.5317, 0.5169, 0.5022, 
0.4431, 0.4283, 0.4136, 0.3988, 0.3840, 
0.3250, 0.3102, 0.2954, 0.2807, 0.2659, 
0.2068, 0.1920, 0.1772, 0.1624, 0.1477, 
0.0886, 0.0738, 0.0591, 0.0443, 0.0343, 
0.0080, 0.0063, 0.0047, 

0.0005, 0.0001, 0.0001, 0.0004, 

0.0044, 0.0045, 0.0177, 0.0310, 0.0443, 
0.1034, 0.1181, 0.1329, 0.1477, 0.1624, 
0.2216, 0.2364, 0.2511, 0.2659, 0.2807, 
0.3397, 0.3545, 0.3693, 0.3840, 0.3988, 
0.4579, 0.4726, 0.4874, 0.5022, 0.5169, 
0.5760, 0.5908, 0.6056, 0.6204, 0.6352, 
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0 . 6499 , 

0 . 7681 , 

0 . 8862 , 

0 . 9855 , 

AIRFY : 

0 . 0093 , 

0 . 0089 , 

0 . 0381 , 

0 . 0546 , 

0 . 0643 , 

0 . 0703 , 

0 . 0707 , 

0 . 0570 , 

0 . 0309 , 

0 . 0033 , 

0 . 0021 , 

0 . 0130 , 

0 . 0468 , 

0 . 0974 , 

0 . 1299 , 

0 . 1411 , 

0 . 1359 , 

0 . 1157 , 

0 . 0864 , 

0 . 0539 , 

0 . 0259 , 

&END 



0 . 6647 , 

0 . 7828 , 

0 . 9010 , 

0 . 9908 , 

0 . 0044 , 

0 . 0135 , 

0 . 0407 , 

0 . 0561 , 

0 . 0652 , 

0 . 0709 , 

0 . 0699 , 

0 . 0540 , 

0 . 0274 , 

0 . 0004 , 

0 . 0033 , 

0 . 0144 , 

0 . 0543 , 

0 . 1025 , 

0 . 1323 , 

0 . 1413 , 

0 . 1340 , 

0 . 1124 , 

0 . 0824 , 

0 . 0498 , 

0 . 0244 , 



0 . 6942 , 

0 . 8124 , 

0 . 9305 , 

0 . 9988 , 

0 . 0000 , 

0 . 0218 , 

0 . 0453 , 

0 . 0589 , 

0 . 0669 , 

0 . 0717 , 

0 . 0675 , 

0 . 0478 , 

0 . 0202 , 

0 . 0001 , 

0 . 0063 , 

0 . 0167 , 

0 . 0680 , 

0 . 1120 , 

0 . 1363 , 

0 . 1411 , 

0 . 1297 , 

0 . 1054 , 

0 . 0744 , 

0 . 0415 , 

0 . 0176 , 



0 . 7090 , 

0 . 8271 , 

0 . 9453 , 

1 . 0000 , 

0 . 0008 , 

0 . 0256 , 

0 . 0474 , 

0 . 0601 , 

0 . 0676 , 

0 . 0719 , 

0 . 0659 , 

0 . 0446 , 

0 . 0165 , 

0 . 0005 , 

0 . 0080 , 

0 . 0168 , 

0 . 0743 , 

0 . 1163 , 

0 . 1379 , 

0 . 1406 , 

0 . 1272 , 

0 . 1018 , 

0 . 0704 , 

0 . 0373 , 

0 . 0124 , 



0 . 7238 , 

0 . 8419 , 

0 . 9601 , 

0 . 9996 , 

0 . 0018 , 

0 . 0291 , 

0 . 0494 , 

0 . 0612 , 

0 . 0683 , 

0 . 0719 , 

0 . 0641 , 

0 . 0413 , 

0 . 0127 , 

0 . 0011 , 

0 . 0097 , 

0 . 0244 , 

0 . 0804 , 

0 . 1204 , 

0 . 1391 , 

0 . 1399 , 

0 . 1246 , 

0 . 0981 , 

0 . 0663 , 

0 . 0331 , 

0 . 0093 , 



0 . 7533 , 

0 . 8714 , 

0 . 9802 , 



0 . 0038 , 

0 . 0353 , 

0 . 0530 , 

0 . 0633 , 

0 . 0697 , 

0 . 0713 , 

0 . 0597 , 

0 . 0344 , 

0 . 0061 , 



0 . 0391 , 

0 . 0920 , 

0 . 1272 , 

0 . 1407 , 

0 . 1375 , 

0 . 1188 , 

0 . 0904 , 

0 . 0580 , 

0 . 0274 , 



0 . 7385 , 

0 . 8567 , 

0 . 9749 , 



0 . 0028 , 

0 . 0323 , 

0 . 0513 , 

0 . 0623 , 

0 . 0690 , 

0 . 0717 , 

0 . 0621 , 

0 . 0379 , 

0 . 0088 , 

0 . 0114 , 

0 . 0319 , 

0 . 0863 , 

0 . 1240 , 

0 . 1401 , 

0 . 1389 , 

0 . 1218 , 

0 . 0943 , 

0 . 0622 , 

0 . 0289 , 



0 . 6795 , 

0 . 7976 , 

0 . 9157 , 

0 . 9956 , 

0 . 0012 , 

0 . 0178 , 

0 . 0431 , 

0 . 0575 , 

0 . 0661 , 

0 . 0714 , 

0 . 0688 , 

0 . 0509 , 

0 . 0238 , 

0 . 0001 , 

0 . 0047 , 

0 . 0157 , 

0 . 0613 , 

0 . 1074 , 

0 . 1345 , 

0 . 1413 , 

0 . 1320 , 

0 . 1090 , 

0 . 0784 , 

0 . 0456 , 

0 . 0219 , 
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APPENDIX F: RVCQ3D CODE INPUT 



A sample namelist input file for RVCQ3D is shown below. Namelist 6 
contains the definition for a linear 6% streamtube contraction which accounts for AVDR 
and the “Quasi-3D” effect in the code. All solutions were run to 6,000 iterations and 
consumed, on the average, only 55 minutes on the NPS CRAY (“Sirius”) and 4 minutes 
on the NAS CRAY (“Von Neumann”). 



'SANGER CONTROLLED DIFFUSION BLADE’ 

&NL1 M=240,N=49JV1TL=60,MIL=108 &end 

&NL2 NSTG=4 JVDT-1 JRS=1 ,EPI=.30,EPJ=.40,CFL=4.0 AV2=0., 

AV4=1. &END 

&NL3 IBCIN=1,IBCEX=1 ,rrMAX=6000aRESTI=04RESTO=l,IRES=MCRNT=10(»00, 
KRM=0 &END 

&NL4AMLE=O.24ALLE=43.03ETE=04»RAT=0.975O > POIN=L0,T0IN=L0,G=L40, 

&END 

&NL5 ILT=2JEDGE=30,RENR=1.0e6JPRNR=.70JPRTR=().9,TW=0.0,VISPWR= 666666, 
CMUTM=14.0 &END 

&NL6 OMEGA=0.0,NBLADE=1 JSfMN=2 &end 
-.3000 3.000 
1.0000 1.000 
1.0000 .9434 
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APPENDIX G: DENSITY RESIDUALS FOR CFD SOLUTIONS 



40 Degree SolulfoM (07o AVDH) 





llerailons 
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APPENDIX H: Cp VERSUS X/C FORTRAN PLOTTING PROGRAM 



cccccccccccccccccccccccccccccccccccccccccccccccc 

C Program To Read The Output From Rvcq3d.F And Grape.F 
C And Generate A Data File For A P/PO Vs Chord Plot 
C Around The Airfoil. Aug.22, 1991 

cccccccccccccccccccccccccccccccccccccccccccccccc 

C IMPLICIT REAL*8(A-H,0-Z) 

DIMENSION Q(300,300, 4),X(300,300),Y(300,300) 

DIMENSION U(300,300),V(300, 300), P(300,300) 

C 

C CALCULA TE THE CRITICAL VELOCITY 

cepe=1.00 
rgas=1.0 
pi=1.0 
ti=1.0 

PRAT=0.685 

rhoi=pi/(rgas*ti) 

ceve=cepe-rgas 

g=1.40 

gP=g + l 

cstar=sqrt(2*g*pi/(gp*rhoi)) 

C 

ISTART=60 

IFINIT=181 

C 

READ(7,*)NI,NJ 
print *,ni,nj 

READ(7,*)((X(I,J),I=1,NI),J=1,NJ), 

& ((Y (I,J),I=1,NI),J=1,NJ) 

READ(7,*)MTL,MIL 

C 

READ(3,*)NI,NJ 

READ(3,*)FSMACH,ALF,RE,TIME 
RE AD(3, *)(((Q(I, J ,K),I= 1 ,NI), J= 1 ,NJ),K= 1,4) 

GAMMA=1.4 
DO 1234 I=1,NI 
DO 4321 J=1,NJ 

U(I,J)=Q(I,J,2)/Q(I,J,1) 

V(I,J)=Q(I,J,3)/Q(I,J,1) 

P(I,J)=(GAMMA-1)*(Q(I,J,4)-.5*Q(I,J,1)*(U(I > J)**2+V(I,J)**2)) 

& *RHOI*CSTAR**2 
4321 CONTINUE 
1234 CONTINUE 
C 

C DETERMINE THE XMIN AND XMAX GRID POINT POSITION 

XMTN=0.0 
XMAX=0.0 

DO 2 I=ISTART,IFINTT 
IF(Xa, 1).LE.XMIN)THEN 
IMIN=I 
PRINT *,IMIN 
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ri n ft 



o o ° r>° onn° o on 



XMIN=X(I, l) 

END IF 

IF(X(I, 1).GE.XMAX)THEN 
EMAX=I 

c PRINT *,IMAX 
XMAX=X(I,1) 

END IF 
2 CONTINUE 

FREE STREAM STA TIC PRESSURE USES THEXM1N GRID POINT 
P0=P(IMIN,NJ)+.5*Q(IMIN ; NJ > 1)*(U(IMIN ( NJ)**2+V(IMIN,NJ)**2) 
& *RHOI*CSTAR**2 

i=imin 

j=nj 

tO=T/rgas*((g-l)*q(ij,4)/q(ij,l)-(g-l)**2/(2*g)* 

& (q(ij,2)**24q(ij,3)**2yq(ij,l)**2) 
t0=t0*cstar**2 

print *,pO,tO,p(inun,nj),q(imin,nj,l),u(imin,nj),v(imin,nj) 

CALCULA TE THE MASS FLOW RA TE FROM THE INPUT DA TA 
NORMALIZED WITH THE INLET AREA 
dmass=pi/(rgas*ti)*(l/prat)**(-l/g)*sqrt(g*(g-l)*rgas*ti/2* 

& ( 1 - 1 /( 1 /prat) * * ((g- 1 )/g))) 

print *,'cal. inlet mass flow rate - ,dmass,'* area' 

dmasse=q(imin,nj,l)*RHOI*u(imin,nj) 
print *,'comp. inlet mass flow rate / area - ,dmasse 

PRINT OUT DOWN STREAM CONDITIONS 
pl=p(l,l) + -5*q(l,l,l)*(u(l,l)**2+v(l,l)**2)*RHOI*CSTAR**2 
tl=l/rgas*((g-l)*q(l,l,4)/q(l,l,l)-(g-l)**2/(2*g)* 

& (q(l,l,2)**2+q(l,l,3)**2)/q(l,l,l)**2) 

tl=tl*CSTAR**2 
c print *, 'downstream condition' 

c print *,pl,tl,p(l,l),q(l,l,l),u(l,l),v(l,l) 

C 

C CALCULATE THE MASS FLOW RATE AT EXIT 

smass=0.0 
ra = 0.0 
ba=0.0 

do 21 j=l,nj-l 

smass=smass+(q( 1 j, 1 )+q( 1 j+ 1 , 1 )) *(u( 1 j )+u( 1 , J+ 1 ) )* 

& RHOI*CSTAR*(abs(y(l j+l)-y(l j)))*0.25 

ra=Ta+(q( 1 j, 1 )+q( 1 j+ 1 , l))*(abs(y( 1 j+ 1 )-y( 1 j)))*0. 5 * 

& RHOI 

ba=ba+abs(y( 1 j+ 1 )-y( 1 j)) 

21 continue 
do 22 j=l,nj-l 

smass=smass+(q(nij,l)+q(nij+l,l))*(u(nij)+u(ni,J+l))* 

& RHOI*CSTAR*(abs(y(ni j+ 1 )-y(ni j)))*0.25 

ra=ra+(q(nl j,l)+q(nij+ 1 , l))*(abs(y(ni j+ 1 )-y(ni j)))*0. 5 
& *RHOI 

ba=ba+abs(y(ni j+ 1 )-y(ni j)) 
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22 continue 



C 

C A VERA GE VELOCITY 

va=smass/ra 
dr=ra/ba 

c print *,Vel. downsteam,vd =',va 

c print *,'ave. density.dr =\dr 

C 

C CALCULA TE THE INPUT INLET VELOCITY 

C NORMALIZED BY CSTAR 

cO=sqrt(g*rgas*ti) 

vin=sqrt(g*(g- 1 ) *rgas* ti/2 *( 1 -( 1/prat)* *(( 1 -g)/g)))/cO 
c print *, 'input inlet velocity / cO - ,vin 
C 

C CALCULA TE THE CORRECTED DOWNSTREAM VELOCITY 
do 23 j=nj,l,-l 
vac=u(nij)*vin/va 
alpha=atan(abs(v(ni j)/u(ni j ))) 
vtc=vac*tan(alpha) 
vc=sqrt(vac**2+vtc**2) 
yy=abs(y(ni,nj)-y( 1 ,nj)) 
c write(24,*)abs(y(nij)-y(ni,nj))/yy,vc 
23 continue 
C 

do 25 j=2,nj 
vac=u(l j)*vinA'a 
alpha=atan(abs(v( 1 j)/u( 1 j))) 
vtc=vac*tan(alpha) 
vc=sqrt(vac * * 2+vtc* * 2) 
c \vrite(24,*)abs(y(lj)-y(ni,nj))/yy,vc 
25 continue 
C 
C 

CHORD=ABS(X(IMAX, 1)-X(IMIN,1)) 

C 

DO 1 I=ISTART,IFINIT 
DIST=ABS(X(1, 1)-X(IMIN, 1 )) 

XS=DIST/CHORD 

CP=1.0+(-P0+P(I,l))/(.5*Q(IMIN,NJ,l)*(U(IMIN,NJ)**2 
& +V(IMIN,NJ)**2)*RHOI*CSTAR**2) 

wcr=sqrt(u(i > l)**2+v(i,l)**2) 

C THE SURFACE POINT Pa, 1 ) ALWAYS EQUAL 0 ? 

PT=P(L1)/P0 

c WRITE(63,*)XS,WCR 
c WRITE(64,*)XS,PT 
WRITE(65,*)XS,CP 
1 CONTINUE 
C 

STOP 

END 



98 



APPENDIX I: UPSTREAM FLOW ANGLE FORTRAN PROGRAM “angle.f” 



ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c 'angler READS FILES 'fort7’ & ’fort.3’ WRITES OUTPUT TO FILE 'fort.10' 
c PROGRAM READS THE OUTPUT FROM RVCQ3D.I AND GRAPE.f AND GENERATES A 
c DATA FILE CONTAINING THE Y-POSITION AND INLET FLOW ANGLE AT A DESIRED 
c LOCATION UPSTREAM OF THE BLADE IN TERMS OF PERCENT CHORD. DATA IS 
c CALCULATED USING SIMPLE LINEAR INTERPOLATION, 
c SEPTEMBER 29, 1994 WRITTEN BY LT. DREW WILLIAMS 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 



real beta(300) 

dimension x(300,300),y(300,300),q(300, 300,4) 
dimension xx(300),yy(300),u(300,300),v(300,300) 
dimension uu(300),w(300) 

open(10) 
rewind(lO) 
pi=3. 14 15926 

c Read in the C-grid created by GRAPE 
read (7,*) ni,nj 

read (7,*) ((x(ij),i=l,ni)j=l,nj), 

& ((y(ij),i=l,ni)j=l,nj) 

read (7,*) mtl.mil 
c write(10,*) x(mtl,l) 

c Read in the Solution file created by RVCQ3D 
read (3,*) ni,nj 

read (3,*) fsmach,alpha,reyn,iters 
read (3,*) (((q(ij,k),i=l,ni)j=l,nj),k=l,4) 

c Calculate the velocities u and v at every grid point 
do 100 i = l,ni 
do 200 j = l,nj 
u(ij)=q(ij,2)/q(ij,l) 
v(ij)=q(ij,3)/q(ij,l) 

200 continue 
100 continue 

c Define the Chord length and then prompt for an upstream location 
chord=x(mtl,l) 

wnte(*,*) 'Input Upstream Chord Location in Percent (i.e., -10, -20,...)' 
c (NOTE: You must be aware of the maximum upstream position available from your 
c GRAPE output. Beyond the maximum upstream location, the values calculated 
c are all the same as that at the maximum.) 
read(*,*) ax 
axperc=ax/100 
vertx=axperc*chord 
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c 

c 



write(10,*) chord 
write/ 10,*) vertx 



c Bracket the vertical x-position 

c (Note: Depending on your grid size, the 7' index may have to be adjusted to 
c ensure the entire inlet flow region is covered.) 
do 300 i = 1,240 
do 400 j = 1 ,nj 
flag=0 

if (x(ij).gt.vertx) jsxl=j 
if (x(ij).lt. vertx) then 
jsx2=j 
flag=2 
goto 500 
endif 
400 continue 
500 continue 

if (flag.gt.l) then 

c Calculate the u and v velocities and compute the inlet angle 
c Linear Interpolation in the form of: ynew = yold + (dy/dx)*dx 
yy(i)=y(i jsx 1 )+(vertx-x(i jsx 1 ))*(y(i jsxl )-y(i jsx2))/ 

& (x(ijsxl)-x(ijsx2)) 

uu/i)=q/ijsxl,2)+/vertx-x/ijsxl))*(q/ijsxl,2)-q/ijsx2,2))/ 

& (x(ijsxl)-x(ijsx2)) 

w(i)=q(ijsxl,3)+/vertx-x(ijsxl))*(q(ijsxl,3)-q(ijsx2,3))/ 

& (x(ijsxl)-x(ijsx2)) 

beta(i)=atan(w(i)/uu(i))* 1 80/pi 
write/ 10,*) yy(i),beta(i) 
endif 

300 continue 
stop 
end 
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